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FOREWORD 


This  document  presents  the  research  performed  under  support  from  the  U.S. 
Air  Force  Office  for  Scientific  Research,  under  Grant  No.  80-0267,  in  the 
period  1980-1984  at  Montana  State  University.  Principal  investigator  was  the 
author.  Dr.  Anthony  Demetriades,  Professor  of  Mechanical  Engineering  at  MSU. 
The  AFOSR  officer  in  charge  was  Captain  Michael  Francis. 

The  objective  of  the  research  was  the  study  of  the  hydrodynamic  stability 
of  a  flat-plate  laminar  boundary  layer  in  the  MSU  Supersonic  Wind  Tunnel  by 
measuring  the  amplification  of  natural  disturbances  in  the  flow.  Initially, 
some  emphasis  was  placed  on  the  additional  effect  of  surface  roughness  on 
stability,  and  a  master's  thesis  by  a  graduate  student  was  prepared  on  the 
flowfield  over  a  rough  wall.  Emphasis  was  gradually  shifted  to  the  smooth 
wall  problem,  which  was  found  more  demanding  and  potentially  more  fruitful 
than  originally  thought. 

The  author  is  indebted  to  AFOSR  for  its  patience  and  support  throughout 
this  work,  and  to  Glenn  McCullough  for  assistance  with  the  model,  the  wind 
tunnel  operation  and  the  gathering  of  the  mean-flow  data. 
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1 .  Introduction  and  Mot ivation 


It  is  common  experience  that  surface  roughness  promotes  transition  to 
turbulence  in  the  laminar  boundary  layer.  Nearly  every  text  on  boundary 
layers  includes  a  summary  of  the  numerous  experiments  done  to  date,  by  which 
transition  to  turbulence  was  found  to  move  upstream  when  the  surface  is 
roughened.  The  body  of  available  literature  is  necessarily  large  because  of 
the  large  variety  of  possible  roughness  geometries  and  their  distribution  on 
the  surface,  i.e.,  of  the  spectrum  of  the  surface  contour.  There  are  several 
prevalent  notions  of  turbulence  generation  by  roughness:  one  is  that  the 
turbulent  wakes  of  a  few  isolated  surface  protrusions  agitate  the  boundary 
layer  into  a  turbulent  state;  another,  that  the  roughness  distorts  the  mean 
flow  field  into  a  hydrodynamically  unstable  shape.  The  latter  view  has 
attracted  attention  especially  when  the  surface  is  uniformly  covered  by 
"distributed"  (statistically  stationary)  roughness  of  height  k  much  smaller 
than  S  (the  layer  thickness). 

As  an  ideal  example  of  a  formal  connection  between  roughness  and 
transition,  one  could  calculate  the  mean  velocity  profile  distortion  due  to 
small-scale,  uniformly  distributed  roughness,  and  then  subject  this  profile  to 
hydrodynamic  stability  analysis;  a  rational  connection  between  the  roughness 
and  transition  would  thus  be  found.  Such  a  task  would  be  arduous  because  of 
the  difficulty  of  the  flowfield  calculation  and  the  need  to  repeat  it  for 
every  conceivable  type  of  roughness.  As  an  alternative,  Reshotko 
(Reference  1)  and  Kendall  (Reference  2)  instead  attempted  to  measure 
experimentally  the  velocity  profile  with  roughness,  with  a  view  of  perhaps 
using  the  measured  profile  as  an  input  to  stability  analysis.  One  could  then 
make  parallel  stability  (e.g.  disturbance  amplification)  measurements,  and 
compare  the  latter  with  stability  characteristics  predicted  from  the  measured 


mean  rough  wall  flowfield. 


The  Reshotko  and  Kendall  tests  were  done  at  low  speeds.  The  work 
described  here  was  planned  as  the  analog  for  supersonic  flows.  Specifically, 
the  purpose  here  was  to  measure  both  the  mean  profile  and  the  amplification 
rates  (stability  diagram)  of  a  supersonic  laminar  boundary  layer,  when  the 
wall  surface  is  rough.  Under  the  best  of  circumstances,  it  was  hoped  that 
eventual  use  of  the  measured  mean  profile  could  be  made  by  stability  theory 
for  amplification  calculations,  and  that  the  stability  characteristics  so 
calculated  would  in  turn  be  compared  with  the  measured  stability 
characteristics.  No  information  exists  to  date  on  the  amplification  of  small 
disturbances  in  a  supersonic  boundary  layer  over  a  rough  wall;  such 
information  would  be  in  any  way  invaluable  toward  the  understanding  of  the 
role  of  roughness  in  promoting  transition.  Thus,  the  data  could  play  a  dual 
role  as  checks  of  the  stability  theory  and  as  practical  guides  to  transition 
prediction. 

At  the  inception  of  the  present  program,  it  was  clearly  understood  that 
previous  knowledge  on  the  supersonic  boundary  layer  stability  with  a  smooth 
wall  should  be  the  necessary  base  on  which  the  measurements  with  roughness 
should  rest.  It  soon  became  apparent  that  such  knowledge  was  overestimated. 
A  survey  of  the  experiments  done  on  smooth  wall  stability  showed  a  number  of 
reports  dealing  with  subsonic  edge  Mach  Number  (e.g.  References  3  and  4),  a 
series  of  experiments  at  1.5  <  <  2.2  (Reference  5)  and  a  rather  heavy 

concentration  at  6  <  Mg  <  8.5  (References  6  through  11  ).  Kendall  (Reference 
12)  made  another  series  of  measurements  at  Mg  ^  3  and  4.5,  but  his 
presentation  deals  mainly  with  the  issue  of  boundary  layer  response  to  the 
free-stream  noise  with  little  information  on  the  disturbance  behavior  within 
the  boundary  layer  especially  at  Mach  3. 
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3  have 


Those  with  some  experience  in  amplification  measurements  at 
given  discouraging  accounts  of  its  suitability  as  a  test-bed  of  linear 
stability  theory.  Laufer  and  Vrebalovich  (Reference  5)  limited  their 
published  account  of  stability  to  ~  1.6  and  2.2  because  "...at  M  >  3  the 
detection  of  self-excited  oscillations  was  much  more  difficult  and  less 
reliable."  Kendall  notes  that  in  his  supersonic  experiments  "fluctuations  of 
all  frequencies  were  observed  to  grow  monotonically  larger  in  the  region  of  a 
boundary  layer  extending  from  the  flat  plate  leading  edge  to  the  predicted 
location  of  instability,  i.e.,  in  a  region  where  no  growth  was  expected" 
(Reference  12,  p.  291).  This  statement  portends  grave  difficulties  for 
stability  experiments  aiming  at  the  observation  of  neutral  boundaries  for 
checking  the  linear  stability  theory  at  or  around  M  -  3.  Such  experiments, 
furthermore,  also  depend  on  amplified  "Tollmien-Schlichting  wave"  observation 
as  a  reliable  indicator  of  ongoing  instability,  and  indeed  the  accidental 
discovery  of  such  waves  by  Schubauer  and  Skramstadt  in  the  1940'$  (Reference 
3)  supplied  the  major  impetus  for  modern-day  stability  research.  Even  at 
hypersonic  speeds,  laminar  instability  waves  are  so  pronounced  that  they  are 
routinely  visible  even  with  unsophisticated  sensors.  This  selectivity  of  the 
boundary  layer  apparently  disappears  at  =  3,  however,  giving  the 
experimenter  no  immediate  evidence  of  disturbance  amplification. 

An  interesting  theoretical  explanation  of  the  exceptional  non-selectivity 
and  low  amplification  in  the  vicinity  of  ■  3  is  supplied  by  Mack  (Reference 
13,  p.  282).  It  turns  out  that  =  3  lies  at  the  minimum  of  curves  one  can 
plot  theoretically  of  maximum  spatial  amplification  rate  vs.  H^.  This  minimum 
marks  the  intersection  of  3-D,  first-mode  amplification  rates,  and  the  rates 
due  to  2-D  second-mode  disturbances.  Thus  ~  3  occupies  a  unique  spot  in 
boundary  layer  stability,  one  which  should  present  difficulties  to  the 


experimentalist  and  the  theoretician  alike. 

As  a  result  of  the  ideas  expressed  above,  the  objective  of  measuring  the 
growth  of  damping  of  natural  disturbances  over  the  smooth  wall,  rather  than 
being  an  initial  tare  measurement,  became  quite  imperative.  The  smooth  wall 
stability  measurement  provided,  in  the  end,  most  of  the  measurements  described 
and  conclusions  reached  here. 

2.  Wind-Tunnel  Facility 

All  measurements  described  here  were  done  in  the  continuous  supersonic 
wind-tunnel  at  HSU  (MSD/SVT)  at  Mach  number  3.0.  A  detailed  description  of 
the  facility  appears  in  Reference  15.  The  relevant  attributes  of  this 
facility  are  its  ability  to  run  for  long  periods  (e.g.  8  hours)  at  constant 
supply  (stagnation)  pressures  and  temperatures,  its  steadiness  and  uniformity 
of  flow,  its  convenience  of  access  to  the  test  section,  its  broad  expanse  of 
optical  view  of  the  flow,  its  automated  probe  control  and  data  acquisition  and 
the  ease  of  controlling  the  sidewall  boundary  layer  transition  zone.  An 
overall  view  of  the  facility  is  shown  in  Figure  1. 

3.  Early  Experiments  with  the  Axi-Svmmetric  Model 

This  program  began  as  a  graduate  student  thesis  experiment  to  look  at  the 
rough  wall  flowfield,  stability  and  transition  on  an  axi-symmetr ic  (ogive- 
cylinder)  model  at  Mach  3.  This  geometry  was  chosen  mainly  to  alleviate 
possible  problems  of  model-wall  interference  common  with  flat-plate  models. 
This  phase  of  the  program  is  presented  in  detail  in  Reference  16. 

The  model  consisted  of  a  20.3  cm  long,  2  cm  diameter  cylinder  attached  to 
an  11.7  cm  long  ogive  with  a  sharp  tip  of  3.2^  half  angle.  The  rear  end  of 
the  ogive  screwed  on  and  blended  smoothly  with  the  front  end  of  the  cylinder 
while  the  latter  was  supported  in  the  back  by  a  sting.  This  model  was  always 
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operated  at  zero  angle  of  attack,  in  the  tunnel  stagnation  pressure  p^  range 
of  200-600  torr  (unit  Reynolds  number  range  20,000  <  Re'  <  60,000/cm),  and 
stagnation  temperature  range  75-125F.  The  model  is  pictured  on  Figure  2. 

The  boundary  layer  flowfield  over  this  model  was  first  examined  when  the 
cylindrical  afterbody  had  a  smooth  surface  (the  ogive  was  always  configured 
with  a  smooth  surface).  The  transition  dependence  on  p^  was  measured,  and  it 
was  next  attempted  to  determine  the  surface  roughness  suitable  for  stability 
measurements  by  introducing  roughness  on  the  cylinder  and  studying  changes  in 
the  transition  location.  Accordingly,  duplicates  of  the  cylindrical  afterbody 
were  built  which  were  covered  by  uniformly  distributed  sand-type  roughness; 
these  roughness  "overlays",  pictured  on  Figure  3,  were  made  of  ordinary  shop 
sandpaper  which  had  earlier  been  subjected  to  a  measurement  of  roughness 
height  by  profilometry  (Reference  17).  The  random  roughness  height  k,  quoted 
below,  was  consequently  known  with  some  accuracy  (as  Reference  17  explains,  k 
represents  the  average  peak-to-valley  height). 

Two  important  findings  emerged.  First,  it  was  extremely  difficult  to 
trip  transition  with  such  sand-type  roughness,  even  though  k  was  gradually 
increased  by  changing  to  coarse  overlays.  For  example,  changes  in  transition 
were  practically  imperceptible  for  60-grit  roughness  (k  “  0.004"  "  0.01  cm). 


For  such  k  we  computed 


«  Weft. 

Re  =  — - —  —  .5^00 


c;  23 


(1) 

(2) 


where  "e"  refers  to  boundary  layer  edge  properties  and  "k"  to  properties  at 
y  =  k.  Feindt  (Reference  18)  claims  that  when  the  former  of  these  two 
Reynolds  numbers  is  about  120,  the  roughness  becomes  effective  as  a  transition 
trip  in  incompressible  flows.  In  view  of  eq.  (1),  Feindt's  criterion  is 
obviously  invalid  at  high  speeds.  This  is  hardly  surprising  for  supersonic 


flows  where  the  density  is  low  near  the  wall  and  where  the  roughness  has  to 
protrude  into  higher  flux  regions  to  be  effective.  The  proper  criterion 
involves  the  Re^^  of  eq.  (2)  as  originally  suggested  by  Schiller  (Reference  19) 
and  Smith  and  Clutter  (Reference  20).  Kendall  (Reference  2)  and  Reshotko 
(Reference  1)  have  re-emphasized  that  Re^^  must  reach  a  value  of  about  100  for 
roughness  to  become  effective.  In  the  present  instance  the  resulting  required 
k  was  computed  to  be  about  0.01"  (0.025  cm). 

Second,  it  became  clear  that  such  large  k  values  would  compound 
difficulties  (already  encountered  for  smaller  k)  of  making  meaningful  boundary 
layer  measurements.  An  average  k  of  0.01”  would  produce  an  average  ratio  of 
k/^  a  1/4,  and  on  occasion  much  larger  for  isolated  large  grit  elements.  But 
unless  k/£  <  <1  the  entire  concept  of  a  uniform,  statistically  homogeneous 

random  surface  collapses  since  k  cannot  be  considered  a  characteristic  length 
constant  from  one  point  on  the  surface  to  another.  Furthermore,  the  solid 
surface  itself  can  no  more  be  defined  when  k  is  so  large.  These  conclusions 
are  obviously  true  for  any  supersonic/hypersonic  flow,  and  therefore  work  with 
the  random-distributed  sand-type  roughness  was  discontinued. 

A  roughness  type  was  next  sought  which  would  eliminate  the  difficulties 
just  described.  The  arrangement  settled  upon,  shown  on  Figure  4,  consisted  of 
a  periodic  pattern  of  alternating  ridges  and  grooves  ("teeth")  machined 
tr ansver sel ly  over  the  entire  length  of  a  cylindrical  afterbody;  the  top  of 
the  ridges  defined  a  cylindrical  surface  coinciding  completely  with  the 
largest  diameter  of  the  ogive  (i.e.  the  grooves  were  recessed  and  the  ridges 
did  not  protrude  from  the  model  surface).  The  ridge  height  was  0.036  cm 
(Figure  4),  consistent  with  the  requirement  on  Re^^  as  per  above.  At  the  same 
time,  however,  the  spacing  between  adjacent  ridges  (=  0.072  cm)  was  made  small 
to  prevent  the  boundary  layer  separating  at  the  top  of  a  ridge  from 
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reattaching  on  the  floor  of  the  adjacent  groove.  It  had  already  been  known 
from  the  work  of  Charwat,  Dewey  and  others  (References  21  and  22)  that  for 
length-to-depth  ratios  of  the  type  shown  on  Figure  4,  the  groove  cavity  is 
"open”,  i.e.  the  layer  separating  from  one  ridge  proceeds  parallel  to  the  flow 
to  reattach  on  the  top  of  the  next.  The  present  choice  of  open-cavity  flow 
would  seem  to  create  a  new  virtual  surface  of  the  model  that  is  defined  by  the 
ridge-tops.  In  turn,  such  a  flow  would  be  free  of  point-to-point  anomalies  in 
the  mean  flow  profiles,  but  at  the  cost  of  also  not  taking  full  advantage  of 
the  ridge  height.  It  will  be  shown  latef  that  the  profile  data  taken  showed 
that  the  flow  over  a  groove  was,  in  fact,  very  similar  to  that  over  a  ridge. 
As  to  the  efficiency  of  this  type  of  roughness  in  tripping  transition,  this 
was  indeed  confirmed  immediately  by  Schlieren  and  pitot  measurements,  as  will 
be  shown  in  Section  4.2. 

At  about  this  time,  it  also  became  clear  that  the  boundary  layer 
development  over  the  ogive-cylinder  had  certain  disadvantages.  For  example, 
the  boundary  layer  growth  was  not  of  the  Blasius  type,  and  especially  at  and 
downstream  of  the  shoulder  the  measured  momentum  Reynolds  numbers  Re  g  were 
too  large  (Figure  5).  Such  behavior  is  typical  of  axi-symmetr ic  ogive- 
cylinder  flows  but  is  not  conducive  to  the  study  of  stability.  This  phase 
having  exhausted  its  usefulness,  the  experiments  were  continued  with  a  flat 
plate  geometry. 

4.  Flowf ield  Measurements 
4.1.  Flat  Plate  Model  Geometry 

Beginning  in  1982,  the  program  was  continued  with  the  design  and 
fabrication  of  a  2-dimensional  sharp-tipped  flat  plate  model  for  the  stability 
measurements.  Like  its  axi-symmetric  predecessor,  this  model,  pictured  on 


Figures  6  and  7,  had  provisions  for  changing  its  top  surface  from  a  smooth  to 
a  roughened  one.  This  was  done  by  the  use  of  interchangeable  inserts,  one  of 
which  was  smooth,  and  the  other  roughened  by  parallel  "teeth"  of  the  same 
geometry  as  used  for  the  axi-symmetric  model  (Figure  6).  All  discussion  will 
henceforth  pertain  to  this  flat  plate  model  only,  which  is  pictured  in  Figures 
6  and  7. 

4.2.  Overall  Flow  Characteristics  and  the  Transition  to  Turbulence 

This  report  deals  with  measurements  done  on  the  flat-plate  model 
described  above  placed  in  the  uniform,  steady  Mach  3  flow  of  the  SWT  test 
section  as  pictured  in  Figure  7.  Before  the  amplification  data  are  discussed, 
however,  it  is  necessary  to  describe  the  overall  features  of  the  flowfield  in 
the  SWT  nozzle  and  the  model  installed  within  it. 

The  overall  flow  features  can  he  seen  in  the  Schlieren  photos  of  Figures 
8  and  9  and  in  the  pitot  profiles  of  the  boundary  layer  included  in  Figures  10 
through  13.  The  Schlieren  picture  shows  clearly  the  boundary  layer  and  also 
the  reflection  of  the  leading-edge  shock  wave.  Transition  to  turbulence  is 
also  shown. 

The  tunnel  noise  environment  is  best  illustrated  via  Figure  13.  The 
boundary  layer  transition  onset  on  the  SWT  interior  sidewalls  has  been  studied 
at  intervals  over  nearly  15  years,  and  it  has  always  depended  on  p^  in  the 
manner  illustrated  on  this  figure.  In  the  present  measurements,  the 
significant  p^  levels  chosen  were  p^  =  350,  475  and  600  torr  (mmHg.).  The 
plate  model  was  installed  so  that  its  leading  edge  lay  13.2  inches  (33.5  cm) 
from  the  nozzle  throat.  According  to  the  graph  of  Figure  15,  this  means  that 
at  Pq  *  350  turbulence  radiation  from  the  sidewalls  never  reached  the  plate 
surface;  on  the  other  hand  at  p^  =  600  the  entire  length  of  the  plate  was 
irradiated  from  the  sidewalls.  The  intermediate  case  of  p_  =  475  torr  is  such 


that  sidewall  radiation  along  Mach  lines  first  reaches  the  plate  at  a  distance 


of  order  3  cm  downstream  of  the  leading  edge. 

Transition  measurements  on  the  plate  were  made  on  several  occasions  in 
this  research  with  the  Schlieren  optics,  pitot  probes  and  the  hot-wire 
anemometer.  In  the  present  instance  transition  onset  was  quantified  using  the 
velocity  profiles  which  are  summarized  on  Figures  8  through  14  (the  flowfield 
implications  of  these  profiles  will  be  explained  in  Section  5.3).  If,  for  the 
moment,  we  concentrate  on  the  agreement  between  the  data  and  the  Blasius 
theory,  it  is  evident  from  Figures  10,  etc.  that  at  a  certain  x  distance  for 
each  of  the  six  case's  shown,  the  data  depart  from  the  Blasius  theory.  This 
"first  departure"  point  (or  x)  was  used  here  to  mark  the  transition  onset. 
These  first-departure  points,  listed  on  Table  I  and  plotted  on  Figure  16,  are 
identified  by  means  of  R,  nominal  Re  and  actual  Re  ;  the  latter  values 

represent  the  momentum  Reynolds  numbers  actually  measured,  and  are  discussed 
in  Section  5.3. 


TABLE  I 

TRANSITION  DATA  SDMMART 


I.  REYNOLDS  NUMBERS 

AT  FIRST  DEPARTURE 

SURFACE 

Pq  (mm) 

Rex  (thousands) 

R«'(Rex)^^^ 

Re«N0M 

Re^ACT 

SMOOTH 

350 

353 

594 

394 

446 

475 

404 

636 

422 

489 

600 

449 

670 

449 

533 

ROUGH 

350 

271 

520 

345 

372 

II.  UPPER  LIMIT  OF  MEASUREMENTS  (FLOW  NOT  YET  FULLY  TURBULENT) 


SMOOTH 

350 

476 

690 

458 

553 

475 

706 

840 

558 

768 

600 

783 

885 

588 

848 

ROUGH 

350 

462 

680 

452 

576 

475 

723 

850 

564 

922 

600 

792 

890 

591 

1032 

The  noteworthy  results  of  Figure  16  are  (a)  transition  as  represented  by 

the  first  departure  is  a  function  of  Re'  for  the  smooth  and  rough  wall  alike, 

(b)  the  ridge-groove  roughness  is  effective  in  moving  transition  forward,  and 

(c)  this  effectiveness  is  pronounced  at  Re'  =  29,000/cm  but  less  so  at 
Re'  ”  56 ,000/ cm.  Note,  for  example,  that  at  Re'  “  29,000/cm  a  decrease  of  15% 
(from  Re^.|.®  400,  say,  to  350)  in  Re^represents  a  30%  decrease  in  the  x  of 
the  transition  location.  Also  to  be  noted  from  Figure  16  is  that  there  is  no 
evidence  from  these  data  that  sidewall  irradiation  makes  any  sudden  changes  in 
the  behavior  of  Re^^  (Re'). 

The  lower  part  of  Table  I  should  dispel  any  illusions  that  the  transition 

process  occurs  at  Reynolds  numbers  as  small  as  indicated  by  the  first 

departure.  The  best  estimates  made  of  C^,  to  be  presented  later,  show  that 
fully  turbulent  flow  was  not  attained  even  at  the  farthest  positions  x 
examined  in  this  experiment.  Especially  if  the  actual  Re^  is  considered,  the 
transition  process  is  still  incomplete  in  the  range  550  <  Re ^  <  1050  and 
depends  on  p^  and  plate  surface  configuration.  The  wetted  distance  x  over 
which  transition  is  in  progress  is  thus  quite  large. 

With  the  preceding  discussion  in  mind,  we  have  attempted  some  comparison 
with  transition  correlations  in  wind-tunnels  compiled  by  Pate  (Reference  23) 


and  more  recently  by  Beckwith  (Reference  24)  and  shown  in  Figures  17  through 
20.  The  comparison  is  somewhat  tenuous  since  the  definition  of  "transition" 
varies  from  one  author  to  the  next,  but  it  is  made  in  order  to  uncover  any 
large  and  fundamental  differences  between  the  MSU/SWT  and  other  facilities. 
There  seems  to  be  no  such  difference.  This  finding  is  important;  it  inspires 
confidence  that  there  is  no  fundamental  reason  why  the  present  transition,  and 
presumably  stability,  behavior  should  be  unique  to  the  present  wind-tunnel. 

TABLE  II 

MATRIX  OF  MEAN-FLOW  MEASUREMENTS 

Mach  number  (nominal):  3.0 

Stagnation  pressures:  600  mm  Hg  abs.  (torr) 

475  mm  Hg  abs . 

350  mm  Hg  abs. 

Stagnation  temperature:  100  F  (560  R) 

(125  F  with  pQ  =  350  mm  only) 

Types  of  surface:  Solid  (smooth  wall) 

Rough  (rough  wall) 

Distance  x  from  L.E.:  0-16  cm  (profiles  taken  at  1  cm  increments) 

Data  Coding  Method:  Data  groups  consist  of  a  four-digit  number  with: 

-First  Digit:  Surface  code:  "3"  for  smooth  surface 

"2"  for  rough  surface 

-Second  Digit:  Pressure  code:  "3"  for  p^  “  bOO  mm 

"2"  for  p^  =  475  mm 
"1"  for  Pp  =  350  mm 
-Third  &  Fourth  Digits:  x(cm) 

(Note:  For  the  stability  measurement  done  with  the  anemometer,  the  third  and 
fourth  digits  refer  to  x  in  tenths  of  an  inch,  e.g.  3347  means  that  the  data 
were  obtained  at  x  =  4.7  inch.) 


Instrumentation  and  Procedures  for  Flowfield  Measurements 
Three  principal  instruments  were  used  for  quantitative  flow  diagnosis  of 
the  mean-average  flowfield.  The  model  surface  temperature  was  measured  by 
contacting  an  iron-constantan  thermocouple  with  the  underside  of  each  insert. 
The  surface  static  pressures  were  measured  with  a  static-pressure  probe 
consisting  of  a  0.022"  dia.  tube  with  0.013"  dia.  holes  drilled  around  its 
circumference  and  located  0.20"  downstream  of  the  closed,  sharpened  end  of  the 
tube.  During  the  measurement  the  tube  was  held  parallel  to  the  flow  and  in 
contact  with  the  surface  and  traversed  slowly  backwards  (see  Figure  21)  while 
the  static  pressure  output  was  recorded  continuously.  Tests  showed  that  the 
spatial  resolution  of  this  probe  was  of  order  5/2;  this  is  satisfactory 
considering  that  the  stream-wise  extent  of  the  mapped  region  was  about  15  cm 
or  100 S ,  and  that  normal  gradients  were  not  expected  on  theoretical  grounds. 

Boundary  layer  profiles  were  recorded  by  traversing  a  0.00A"-dia.  pitot 
probe  normal  to  the  surface.  Based  on  an  average  S  of  0.15  cm,  the  spatial 
resolution  of  this  probe  was  of  order  1:15.  The  SWT  probe  actuator  system 
served  both  to  suspend  the  probe  and  to  move  it  in  an  advance-pause-measure 
sequence.  A  schematic  of  the  system  is  shown  on  Figure  22.  The  pitot  signal 
was  translated  from  pressure  to  voltage  by  a  Kulite  0-5  psia  pressure 
transducer,  so  that  at  each  point  of  measurement  two  voltages  were  auto¬ 
matically  recorded,  one  corresponding  to  distance  above  the  surface.  These 
voltages  were  converted  to  digits  by  a  Spectral  Dynamics  SD133  A/D  Converter 
and  stored  on  cassette  by  a  Texas  Instruments  Silent  700  computer  terminal. 

The  mean-flow  m.atrix  of  measurements  is  shown  on  Table  II.  For  each  of 
the  six  combinations  of  3  stagnation  pressures  and  2  types  of  surface,  the 
plate  surface  pressure  and  temperature  was  determined  and  boundary  layer 
profiles  spaced  1  cm  apart  were  obtained  with  the  pitot  probe,  each  profile 


I 
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typically  containing  60  point  measurements-  Several  profiles  were  taken  in 
the  transitional  zone  as  well.  All  data  were  taken  at  the  mid-span  position 
of  the  plate.  Each  profile  was  coded  with  a  four-digit  number  as  explained  on 
Table  II. 

During  the  measurements,  it  was  discovered  that  the  boundary  layer  growth 
became  anomalous  at  x  >  6  cm  when  p^  was  475  and  350  mm.  Since  the  SWT 
sidewall  was  laminar  at  these  pressures  at  the  point  where  it  intersected  the 
plate  leading  edge,  it  was  determined  that  the  latter  caused  sidewall  boundary 
layer  separation  and  set  up  a  system  of  waves  converging  at  x  =  6  cm  on  the 
top  plate  surface.  The  problem  was  corrected  by  limited  tripping  of  the 
sidewall  boundary  layer  ahead  of  its  intersection  with  the  plate  leading  edge. 
Details  of  this  scheme  are  given  in  Appendix  A. 

The  mean-flow  data  were  reduced  by  the  standard  technique  of  combining 
the  pitot,  static  pressure  and  wall  temperature  measurements.  At  each 
position  X  chosen,  the  static  pressure  was  known  from  the  continuous  static- 
probe  traverses  described  above.  The  total  temperature  had  been  measured  at 
only  one  point,  but  this  was  considered  adequate  because  of  the  expected 
uniformity  of  T^  along  the  plate  and  the  known  insensitivity  of  the  data 
reduction  results  on  T^.  The  measured  T^  was  assumed  valid  at  each  point  on 
the  surface  and  each  p^.  The  distribution  of  T^  through  the  layer  was  assumed 
linear;  T^^  at  the  layer  edge  was  found  to  be  equal  to  T^. 

The  data  reduction  programs  did  not  account  for  the  lietortion  of  the 
pitot  probe  reading  due  to  viscosity  (Reference  25).  The  viscous  effect 
should  be  most  pronounced  at  the  surface,  where  the  minimum  diameter-based 
pitot  probe  Reynolds  number  was  typically  about  10.  A  viscous  error  of 
several  percent  is  possible  in  this  case  for  the  data  taken  right  next  to  the 
surface.  On  the  other  hand,  such  data  should  be  in  any  way  rather  unreliable 
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because  of  probe-surface  interference,  and  their  dubious  validity  did  not 
justify  viscous  correction  (see  below  for  the  handling  of  the  first  few  points 
next  to  the  wall). 

The  boundary  layer  thickness  was  determined  by  inspection  of  preliminary 
plots  of  each  profile.  As  always,  the  choice  of  S  is  somewhat  subjective  in 
boundaries  of  laminar  flow  which  are  diffuse,  a  difficulty  compounded  by  the 
finite  probe  size.  Generally  the  method  of  Kendall  (Reference  26)  was 
followed,  and  onceS  was  found  then  final  data  reduction  of  the  profile  was 
made. 

The  local  "edge"  properties  (subscript  "e")  were  found  from  the  last 
profile  point,  which  usually  lay  at  y  —  l.sS  .  The  flow  within  0  <  y  <  1.35 
was  invariably  uniform  and  constant,  so  that  this  method  eliminated  the 
propagation  of  errors  in  finding  8,  to  the  determination  of  the  edge 
properties. 

During  the  measurement,  the  pitot  probe  tip  was  always  observed  with  a 
microscope  to  ensure  that  the  first  recorded  datum  of  each  profile  was  taken 
with  the  probe  touching  the  surface,  and  the  second  with  the  probe  off  the 
surface.  Because  of  the  miniscule  motions  involved  and  the  finite  precision 
of  the  hardware,  more  points  than  one  were  always  recorded  with  the  probe 
touching  the  surface.  This  could  be  confirmed  with  the  microscope,  however, 
and  a  tag  notifying  the  computer  of  such  "bad"  points  was  included  in  the  data 
analysis.  For  the  rough  wall  measurements,  the  "surface"  terminology  must  be 
replaced  with  "top  of  ridge";  note  Figure  6  in  this  connection. 

The  data  thus  reduced  for  the  laminar  portion  of  the  boundary  layer 
produced  velocity  profiles  which  generally  agreed  well  with  the  Blasius 
theory,  as  will  be  shortly  demonstrated,  except  for  a  small  displacement  along 
the  y  axis.  This  displacement  ranged  from  0  to  10  mils  in  extreme  cases,  and 


remains  unexplained  at  this  writing-  To  account  for  this  problem,  the  data 
were  first  plotted  to  measure  this  displacement  ("offset")  for  each  profile, 
and  then  reduced  for  a  second  time  by  adding  the  appropriate  offset  to  the 
measured  vertical  distance  y.  This  came  to  be  known  as  the  "offset"  method  of 
data  reduction.  The  offset  method  had  the  added  feature  that  anomalous  points 
near  the  surface  were  eliminated  while  computing  the  momentum  thickness  and 
also  during  final  plotting  of  the  profiles. 

A  special  problem  arising  in  the  rough  wall  measurement  concerned  the 
definition  of  the  "surface".  As  already  mentioned,  the  succession  of  two- 
dimensional  cavities  formed  by  the  ridge-groove  roughness  was  ignored  by  the 
flow,  in  the  main,  because  each  cavity  was  closed.  Pitot  profiles  taken 
directly  above  the  ridges  gave  velocity  profiles  identical  to  those  taken  over 
the  grooves  (e.g.  midway  between  two  successive  ridges).  Thus,  the  plane 
formed  by  the  tops  of  the  ridges  formed  the  "surface"  in  the  rough  wall 
measurements  (see  Figure  23). 

The  management  of  the  data  is  outlined  in  Appendix  B.  The  main  inputs 
consisted  of  tabulations  of  pitot  readings  vs.  y  into  files  called  PITOXXXX  or 
PITXXXX  where  XXXX  is  the  code  given  in  Table  II.  The  principal  reduction 
program  was  LAMBL2  which  received  the  PITOXXXX  files  and  produced  printouts  of 
the  results  and,  if  requested,  prepared  and  stored  files  GRAFXXXX  suitable  for 
graphics  plots  of  these  results.  Summary  graphics  plots  of  integral 
properties  vs.  x  were  also  possible  by  using  files  called  DOTALXX,  where  here 
XX  refers  to  the  first  two  digits  of  the  code  of  Table  II. 

5.  Results  of  the  Flowf ield  Measurements 
5.1.  Plate  Surface  Temperature 

The  surface  temperature  measurement  gave  a  recovery  factor  of  0.937 
0.007  for  all  conditions  of  p^  and  surface  configurations. 


5.2.  Plate  Surface  Pressure 

The  surface  static  pressure  distribution  along  the  plate  centerline  is 
shown  on  Figure  21.  There  are  discernible  small  pressure  gradients  toward  the 
trailing  edge  where  stability  data  were  not  taken  in  anyway;  otherwise  the 
pressure  is  nearly  constant,  insensitive  to  the  surface  roughness  and  such  as 
expected  from  isentropic  expansion  calculations. 

It  is  of  some  interest  to  examine  the  surface  pressure  distribution  over 
the  rough  wall  in  the  greatly  magnified  detail  of  Figure  24.  The  trace  shown 
has  been  obtained  during  the  static  probe  traverse  (see  Section  4.3),  has  been 
copied  directly  from  the  x-y  plotter  record,  and  is  a  magnified  portion  of  the 
static  pressure  distribution  in  the  neighborhood  of  x  =  12  -  13  cm  for  the 
rough  wall  at  p^  =  600  torr.  For  this  case  the  boundary  layer  included  in  the 
graph  is  turbulent;  note  the  relative  dimensions  mentioned  on  the  figure  (and 
refer  to  Figure  8  for  a  photographic  view).  The  remarkable  feature  of  this 
trace  is  its  periodicity;  there  appear  two  prominent  wavelengths,  the  larger 
(0.32  cm)  being  of  order  S  ,  the  smaller  one  (0.08  cm)  being  very  close  to  the 
spacing  between  ridges,  i.e.  to  the  roughness  wavelength.  The  smaller 
wavelength  is  undoubtedly  associated  with  the  roughness,  i.e.  with  an 
equivalent  wavy-wall  effect.  The  origin  of  the  larger  wave  is  not  known.  In 
either  case  the  periodicity  in  the  pressure  is  rather  small,  of  order  0.62  for 
the  larger,  and  an  estimated  0.052  for  the  smaller  wave.  It  should  be 
stressed  that  the  pattern  of  Figure  24  disappeared  when  the  boundary  layer  was 
laminar,  in  agreement  with  the  notion  that  pressure  features  are  accented  when 
the  flow  over  the  generating  surface  features  increases  in  speed.  No  account 
of  this  waviness  was  therefore  taken  in  the  analysis  of  the  stability  data. 


5.3.  Boundary  Laver  Development 

The  profiles  of  flow  properties  through  the  boundary  layer  are  shown  on 


Figures  10  through  14;  integral  properties  are  shown  on  Figures  25  through  27 
and  tabulated  in  Table  III. 


These  data,  together  with  the  preview  afforded  by  Figures  8,  9  etc. 
collectively  support  all  preconceived  ideas  about  a  normal  laminar  boundary 
layer  ending  into  a  transitional  zone.  For  the  smooth  wall  (e.g.  Figure  13) 
the  velocity  profile  is  in  agreement  with  the  Blasius  theory  (Reference  27)  up 
to  the  "first  departure"  point  already  discussed  in  Section  4.2.  Beyond  that 
point  an  increasing  disparity  appears  between  the  theory  and  the  data  as  x 
increases.  The  agreement  between  the  theory  and  the  laminar  data  is  most 
convincing;  a  few  isolated  exceptions  occur  near  the  leading  edge  where  the 
laminar  boundary  layer  is  so  thin  that  some  probe-flow  interference  is 
expected.  All  profiles  show  no  data  very  near  the  wall  (y/5<  0.2)  for  the 
same  reason.  Note  that  the  data  are  shown  in  terms  of  the  transformed 
coordinate  y  normalized  with  the  measured  momentum  thickness: 


where  all  quantities  in  the  integrands  were  drawn  from  the  measurements. 

The  integral  properties,  etc.  from  Table  III  have  been  compared  with 


their  theoretical  counterparts: 

5  (i^)  ‘  Ma)  (4) 

Ue  “ 

e  =  0.644  ^  (5) 

=  Re'5^  =  •o4X  (6) 

Re  0  =  Re'e’’  =  o.  441  X  (7) 

Ree  -  0,464  =  o.  664R.  (8) 

where  for  the  three  latter  formulas  we  have  used  y  ®  1.4,  ■  3  and  C  *  1. 

Eq.  (8)  has  been  plotted  on  Figure  13  and  14,  and  it  is  seen  that  there  is  a 


TABLE  III.  FLOMFIELD  INTEGRAL  PROPERTIES 
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systematic  and  growing  departure  between  this  formula  and  the  data  as  x 
increases,  which  is  apparently  independent  of  p^^  (and  Re').  Since  Reg  will  be 
needed  later  for  the  description  of  the  stability  behavior,  its  departure  from 
eq.  (8)  was  formalized  by  a  least-squares  fit  of  the  points  shown  on 


Figure  13: 
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where  "ACT"  and  "NOM"  refer  to  actual  (measured)  and  nominal  respectively,  and 

where  the  nominal  Re^  is  given  by  eq.  (8).  The  coefficients  aj,  the 

smooth  wall  were  aj  =  -6.609E-3  and  a2  “  4J15E-4. 

For  the  rough  wall,  the  rather  surprising  fact  emerges  that  the  velocity 

profile  shows  no  characteristic  "signature"  of  the  roughness.  This  can  be 

verified  from  Figures  such  as  14,  etc.  where  it  is  seen  that  any  differences 

between  the  rough  wall  profiles  and  the  Blasius  theory  or  the  smooth  wall  data 

were  very  small  and  irregular  so  long  as  the  boundary  layer  remained  laminar. 

Much  effort  was  devoted  unsuccessfully  in  plotting  various  aspects  of  the  data 

in  order  to  find  any  systematic  differences  between  rough  and  smooth  walls. 

The  best  that  could  be  done  is  pictured  on  Figure  28.  Here  the  abscissa  is 

again  in  terms  of  ~  ,  and  its  range  covers  fully  the  boundary  layer  width 
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As  this  Figure  shows,  for  low-speed  flows  with  a  sandpaper  surface  Kendall 
found  a  positive  increment  near  the  wall,  meaning  that  the  roughness  acted  to 
decelerate  the  flow.  The  data  points  from  this  experiment  indicate  a  certain 
acceleration.  However,  one  should  note  the  data  scatter;  furthermore,  the 
data  shown  are  the  end  products  of  a  heavy  editing  process  (other  data  were 


much  too  scattered  for  inclusion)  and  are  shown  as  a  rather  biased  attempt  at 
comparison.  Mo  conclusions  are  thus  recommended  on  the  basis  of  this  figure, 
although  it  is  clear  that  the  present  data  lie  sufficiently  close  to  the  wall 
to  divulge  any  chance  distortion  there  due  to  the  roughness.* 

The  momentum  Reynolds  number  for  the  rough  wall,  according  to  Figure  14, 
departs  from  the  smooth  wall  theory,  as  already  observed  for  the  smooth  wall 
in  Figure  13  and  eq.  (9).  The  actual  (measured)  Re  for  the  rough  wall  was 
therefore  fitted  with  a  formula  like  that  of  eq.  (9),  except  that  in  the 
present  case  aj  =  -0.1354,  ^2  “  8.741E-4. 

6.  Stability  Measurements 

6.1.  Instrumentation  and  Data  Acouisitioo  Procedure 

The  stability  data  were  extracted  from  the  AC  signal  of  a  constant- 
current  hot-wire  anemometer  probe  such  as  pictured  on  Figure  29.  By 
exercising  extreme  care  in  using  these  probes,  such  as  shielding  them  during 
tunnel  start  or  stop,  only  one  probe  (No.  9)  failed  during  the  many  8-hour 
days  necessary  to  accumulate  the  data,  so  that  a  total  of  only  two  probes 
sufficed  for  the  measurements.  The  geometrical  and  operating  characteristics 
of  the  probes  are  shown  on  Table  IV. 

*The  curve  due  to  Kendall  shown  on  Figure  28  was  deduced  by  him  from  a  single 
boundary  layer  profile  with  and  without  roughness.  The  scatter  of  the  present 
data,  on  the  other  hand,  arises  from  the  superposition  of  several  profiles; 
there  is  little  scatter  within  any  given  profile.  The  issue  therefore  is 
whether  Kendall's  increment  is  self-similar  and  universal  at  low  speeds,  for 
the  type  of  roughness  he  employed. 


TABLE  IV 

HOT-WIRE  ANEMOMETER  CHARACTERISTICS 


PROBE  NO. 

5 

9 

Diameter,  microinch 

20 

20 

Material 

Ft  lOZ  Rh 

Pt  lOZ  Rh 

Resistance  at  0  C  (ohms) 

78.57 

70.8 

Aspect  Ratio 

175 

158 

Heating  Current,  ma 

3.4 

3.03 

Overheat,  percent 

50 

49 

Amplifier  gain 

50 

75 

Time  constant  (microsec.) 

30 

30 

High  cut-off  (khz) 

500 

1000 

Low  cut-off  (HZ) 

3 

3 

The  function  of  the  measurement  was  to  make  a  permanent  analog  record  of 
the  probe  output  at  each  of  a  large  number  of  points  x,  at  constant  y/<S  ,  for 
each  combination  of  p^  and  surface  condition  (rough  or  smooth).  Such  a 
permanent  record  mas  made  possible  with  the  equipment  diagramed  on  Figure  30. 
The  AC  probe  output  was  recorded  on  a  wideband-FM  channel  of  a  seven  channel 
analog  recorder  (Honeywell  7600)  at  120  ips  (Channel  7),  which  had  a  response 
to  400  KHZ;  a  direct  channel  (Channel  2)  with  a  response  to  2  MHZ  also 
received  this  signal,  for  subsequent  study  of  high-frequency  (0.4-2  MHZ) 
phenomena  of  possible  interest.  The  signal  was  monitored  in  real  time  by  the 
devices  shown  on  the  left  on  Figure  30,  while  in  the  upper  right  of  the  Figure 
is  shown  equipment  by  which  the  probe  position  was  controlled,  monitored  and 
recorded.  Channel  4  of  the  recorder  was  used  for  a  20  KHZ  tone  (see  lower 
left  hand  side  of  Figure  30)  which,  on  playback,  was  used  as  a  "valid  data  on" 
command  for  data  reduction.  Finally,  equipment  shown  on  the  lower  right  hand 
side  of  the  figure  was  used  to  ensure  that  the  recording  process  went 
smoothly,  in  real  time. 
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The  measurement  consisted  of  first  positioning  the  probe  tip  at  the 
desired  x  and  y/6  (the  choice  of  the  latter  will  be  explained  in  Gection  6.2) 
and  recording  the  probe  signal  for  10  seconds;  then  the  probe  was  moved  to  a 
new  X,  keeping  y/5  the  same,  and  the  process  repeated.  Each  10-second  signal 
burst  was  called  a  "data  group"  and  logged  as  per  the  rules  of  Table  II.  The 
collection  of  groups  for  each  particular  p^  and  surface  condition  is  called  a 
"data  set".  Data  sets  were  taken  not  only  inside  but  also  outside  the 
boundary  layer,  the  latter  by  holding  the  probe  0.25"  (0.625  cm)  above  the 
plate  surface.  The  stations  x  were  separated  by  0.1"  (0.25  cm)  and  ranged 
from  about  0.2"  downstream  of  the  L.E.  to  a  point  well  into  the  transitional 
regime . 

The  data  are  summarized  in  Table  V.  The  "BL"  designation  refers  to  the 
data  taken  within  the  boundary  layer,  while  "FS"  refers  to  those  taken  in  the 
free  stream,  as  just  explained  above;  there  are  thus  2  x  3  x  2  =  12  sets  of 
data  altogether.  The  number  of  x  positions,  i.e.  data  groups,  is  seen  to  vary 
from  one  set  to  another.  For  example,  for  p^  =  350  torr  with  the  rough  wall, 
there  were  12  +  17  =  29  groups  recorded  in  the  free  stream  outside  the  layer 
in  the  "forward"  probe  position  with  probe  No.  5,  and  another  14  +  22  =  26 
groups  recorded  in  the  "back"  position  with  probe  No.  9.  Thus  there  was  a 
plate  lenth  equal  to  (29  +  26)  /  10  =  5.5  inches  of  flow  covered  at  0.1" 
intervals.  In  all,  690  data  groups  were  recorded  on  20  reels  totaling  about 
210,000  feet  of  tape. 

It  is  important  to  note  that  the  hot-wire  anemometer  responds  jointly  to 
fluctuations  in  the  fluid  speed,  its  temperature,  density  and  pressure.  The 
process  by  which  the  latter  fluctuations  are  extracted  from  the  wire  AC 
voltage  is  called  "modal  analysis".  In  practice  (References  5,  8  and  12  for 
example)  modal  analysis  is  put  aside  in  stability  experiments  because  of  its 
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great  complexity,  because  of  the  theoretically-confirmed  insensitivity  of  the 
stability  to  the  precise  mode  of  fluctuation  (Reference  13)  and  because  of 
recent  experimental  confirmation  of  such  insensitivity  by  Stetson  et  al 
(Reference  10).  Therefore,  in  this  work  the  quantity  A(f;x),  while  in  reality 
the  spectral  density  of  the  AC  anemometer  voltage  output,  is  equated  to  the 
rms  spectral  density  of  a  typical  fluctuation.  Spectral  densities  or 
amplitudes  were  in  anyway  unimportant  in  this  work,  and  only  the  amplitude 
change  was  studied.  The  complex  and  lengthy  process  of  computing  this  change 
(i.e.  the  amplification  rates)  from  the  anemometer  signal  is  explained  in 
Appendices  C  and  D. 

6.2.  Initial  Observations 

Initial  observations  with  the  hot-wire  anemometer  confirmed  that  the  zone 
of  transition  to  turbulence  varied  in  the  manner  shown  on  Figure  16.  In  the 
free  stream  the  rms  output  of  the  wire  was  generally  so  low  as  to  be  almost 
indistinguishable  from  the  electronic  noise.  The  wire  signal  always  increased 
in  the  laminar  boundary  layer,  however.  Traces  of  the  wideband  signal  vs.  y/5 
are  shown  on  Figure  31.  In  the  range  0.2  <  y/S  <  0.9  this  increase  consists 
of  a  "double-peaked"  curve  which  appeared  to  maintain  a  self-similar  form 
along  x,  independent  of  p^  and  surface  condition.  Of  these  two  peaks,  the  one 
occurring  at  about  y/g  =  0.45  was  slightly  larger  than  the  other,  occuring  at 
about  y/5  =  0.68.  Similarity  of  this  profile  stopped  at  the  first  departure 
point . 

At  low  speeds,  Schlichting's  theory  (Reference  28)  and  wind-tunnel  data 
(Reference  3)  suggest  such  a  double-peaked  distribution  because  of  shifting 
phase  of  neutral  oscillations  across  the  boundary  layer.  The  experimental 
record  at  hypersonic  speeds,  on  the  other  extreme,  speaks  of  a  similarly 
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"noisy"  laminar  boundary  layer  profile  but  with  a  single  peak  (References  6, 
etc.).  A  double  peak  is  noticed  in  the  early  data  shown  at  M  =  2.2  by  JPL 
(Reference  5),  but  more  recent  JPL  measurements  by  Kendall  in  the  2  <  <  5 

range  indicate  a  double-peak  only  at  =  3.  It  seems,  therefore,  that  in  the 
supersonic/hypersonic  regime  this  double-peaked  wideband  fluctuation  profile 
is  specific  to  or  around  Mach  3,  although  for  reasons  presently  unknown. 

Another  peculiarity  obvious  from  Figure  31  is  that  the  wideband  peak 
signals  do  not  always  intensify  along  x,  and  that  these  signals  are  already 
large  at  very  small  x,  of  order  1  cm.  One  possible  answer  is  that 
considerable  amplification  occurs  very  near  the  leading  edge  (say  at  x  <  1  cm) 
where  linear  stability  theory,  however,  predicts  nothing  but  damping  (Re^  at 
distances  so  close  to  the  leading  edge  are  of  order  100).  Coupled  with  the 
greatly  suppressed  fluctuation  levels  in  the  free  stream,  as  mentioned  above, 
this  implies  a  sudden  large  increase  of  the  disturbance  level  as  the 
fluctuations  first  enter  the  layer  very  near  the  leading  edge.  We  shall 
return  to  this  subject  later  from  another  direction,  when  we  discuss  the  total 
amplifications  measured  in  the  boundary  layer. 

All  measurements  reported  hereafter  were  taken  at  y/5  =  0.68,  that  is  on 
the  energy  peak  of  Figure  31  which  was  farther  from  the  wall.  Using  this  peak 
as  a  landmark  simplified  the  task  of  keeping  ylh  constant  from  one  x  to 
another.  Use  of  the  peak  at  y/S  =  0.45  was  rejected  for  fear  of  increased 
risks  of  wire  breakage  and  wall-probe  interference. 

The  data  taken  outside  the  boundary  layer  (sets  51,  52...,  61,  — )  were 
used  to  scrutinize  the  effect  of  Pq  on  the  stream  turbulence,  via  Figure  32. 
Even  though  the  signal  was  not  modally  analyzed,  it  is  clear  from  this  figure 
that  there  exist  two  general  levels  of  stream  turbulence  intensity,  both  of 
which  incease  with  Re^^.  At  the  top  level,  marked  "with  B.L.  radiation"  in  the 


figure  one  finds  the  data  recorded  when  the  wire  was  immersed  in  the  sound 
radiation  field  from  turbulent  sidewall  boundary  layers.  The  lower  level, 
marked  "w/o  B.L.  radiation"  corresponds  to  laminar  sidewall  boundary  layers. 
The  figure  confirms  an  earlier  projection,  based  on  Figure  15,  about  the 
location  of  sidewall  transition  (see  Section  4.2).  These  "free  stream" 
wideband  signals  range  from  0.015  to  0.04  vrms  compared  to  the  boundary  layer 
signals  (Figures  33  and  34)  which  range  from  about  0.06  to  1.6  vrms.  The 
latter  signals  also  scale,  approximately,  with  Re^^  as  shown  on  Figures  35  and 
36.  This  is  a  most  important  result,  since  it  says  that  the  sidewall 
radiation  has  little,  if  any,  systematic  effect  on  the  growth  of  disturbance 
level  in  the  boundary  layer. 

In  all  circumstances  the  wideband  rms  signal  along  the  boundary  layer  was 
found  to  increase  considerally  some  distance  before  the  point  of  first 
departure.  Figures  37  to  42  have  been  prepared  to  set  this  increase  in 
context  with  observed  profile  shapes  and  friction  coefficients  discussed  in 
Section  5.  It  is  seen  that  this  increase  amounts  to  a  factor  of  3  in  typical 
cases,  and  it  starts  much  closer  to  the  leading  edge  than  sometimes  thought 
especially  if  one  keeps  in  mind  the  relation  x  Re^  .  For  example,  for  the 
rough  wall  at  p^  =  350  torr  (Figure  40)  the  x  at  which  the  wideband  rms  starts 
increasing  is  half  the  x  at  the  first  departure.  Wideband  rms  observations 
alone  would  thus  be  most  unsuitable  for  locating  the  transition  zone. 

6.3.  Results  of  the  Stability  Measurements 
6.3.1  Spectra  of  the  Fluctuations 

The  spectra  of  the  fluctuations  are  shown  typically  in  Figures  43  and  44, 
to  an  upper  limit  of  120  KHZ.  In  accord  with  expectations  from  linear 
stability,  the  spectral  density  seems  to  decrease,  generally,  at  the  high 
frequencies.  The  rough  wall  data  appear,  too,  to  be  much  more  "selective" 
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than  the  smooth  wall  data,  with  peaks  in  the  10-40  KHZ  range.  In  contrast  to 
results  presented  at  subsonic,  low  supersonic  and  hypersonic  speeds  in 
References  3,  5,  6,  etc.  the  present  smooth  wall  results  are  quite 
unselect ive.  This  is  so  despite  the  practice  in  most  earlier  works  to  present 
the  mean-square  spectral  density  A^(f),  as  opposed  to  the  rms  density  A(f) 
shown  on  Figures  43  and  44  which,  incidentally,  is  therein  labeled  "amplifier 
output."  Figure  45  shows  a  comparison  between  the  spectra  in  the  boundary 
layer  and  those  in  the  free  stream  under  the  same  conditions. 

A  three-dimensional  view  of  spectrum  development  for  the  six  cases  under 
study  is  shown  on  Figure  46.  The  selectivity  of  the  rough  wall  boundary  layer 
becomes  very  obvious  in  these  plots.  Note  that  the  plots  extend  into  the 
boundary  layer  transition  zone. 

As  mentioned  above,  the  smooth  wall  data  had  very  low  selectivity.  Even 
so,  peaks  in  their  spectra  could  be  faintly  identified  toward  the  end  of  the 
laminar  range.  Figure  47  shows  the  position,  on  the  stability  diagram  (F, 
Re^),  of  all  the  spectrum  peaks  thus  observed.  The  locus  is  often  called 
"maximum  amplification"  line  in  the  literature,  which  should  not  be  confused 
with  the  "maximum  amplification  rate"  line  in  the  (  F,  Re^  )  space.  The  data 
shown  on  this  figure  were  taken  directly  from  the  spectrum  peaks  found  from 
figures  such  as  45,  etc. 

The  point  made  by  Figure  47  is  that,  first,  the  maximum  amplitude  line  is 
quite  independent  of  Re';  second,  there  seems  to  be  no  effect  of  the  rough¬ 
ness;  third,  the  data  agree  with  those  of  Laufer  (Reference  5)  in  that  they 
form  with  the  latter  a  logical  progression  in  the  range  Mg  =  1.5  -  3.  In  this 
respect,  note  that  the  agreement  improves  when  the  actual  Re^  is  considered. 


6.3.2.  Characteristic  Wavelengths 

At  this  point  it  is  appropriate  to  remark  about  an  interesting 
generalization  implied  by  Figure  Motivated  by  similarities  between  the 
incompressible  and  compressible  boundary  layers  when  coordinates  in  the  latter 
are  compressed  by  the  compressibility  transformation,  Laufer  and  Vrebalovich 
found  (Reference  5,  Fig.  32)  that  their  neutral  amplification  data  coincided 
with  those  at  M  =  0,  by  Schubauer,  if  F  is  plotted  vs. 


(11) 


A  coincidence  similar  to  that  indicated  in  Reference  5  for  the  neutral 
branches  is  bound  to  occur  also  for  the  maximum  amplitude  line  which  follows 
closely  the  upper  neutral  branch,  inviting  an  inclusion  of  the  Schubauer- 
Skramstad  data  in  Figure  47.  We  can,  however,  do  better.  The  upper  neutral 
branch  for  the  amplified  region  observed  in  References  3  and  3  is  very  nearly 


a  hyperbola,  so  that  in  the  correlation  scheme  of  Reference  5  the  maximum 
amplitude  should  occur  along  a  line  of  the  type 
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which  in  turn  gives 

i  (13) 

where  A  is  the  disturbance  wavelength  based  on  the  edge  fluid  (not  the  phase) 


velocity.  Thus  the  Laufer-Vrebalovich  correlation  implies  that  the  wavelength 
is  independent  of  Mg,  a  point  which  can  now  be  formally  checked  by  also 


including  the  present  data.  In  fact  one  can  also  make  a  comparison  with  the 


hypersonic  data,  where  has  been  found  to  be  fairly  independent  of  Re  as 

s  s 

well. 


This  comparison,  where  A  vas  computed  from  the  edge  velocity  and  the 
measured  peak  frequencies,  is  shown  on  Figure  48  for  all  stability  data.  A 


few  points  obtained  by  Kendall  (Reference  12)  at  «  7.7  are  not  shown,  but 
these  data  fall  almost  exactly  on  top  of  the  other  hypersonic  results 
included.  Belov  the  bottom  margin  of  the  figure  lie  "higher  harmonic"  results 
about  which  there  is  good  agreement  among  References  6,  9,  and  12.  In  passing 
it  must  also  be  noted  that  (a)  the  present  data  are  plotted  using  the  actual, 
measured  momentum  thickness  (see  previous  discussion),  (b)both  smooth  wall  and 
rough  wall  data  from  the  present  test  are  included  and  (c)  the  M  *  7  data  of 
this  author  represent  three  different  experiments  and  include  hot-  and  cold- 

wall  results  (to  T„/T-  =  0.4). 

V  o 

It  is  evident  from  this  figure  that  the  prominent  wavelengths  for 
0  <  Mg  <  3  are  of  order  30  ,  while  at  hypersonic  speeds  they  are  much  smaller 
in  terms  of  S  (  ^'=>22,  as  previously  found  in  References  6 , 9  etc.).  The 
0  <  Mg  <  3  results  seem  unaffected  by  roughness,  while  the  M  7  results  seem 
unaffected  by  moderate  cooling.  Thus  this  large  difference  foretells  of 
fundamentally  different  instability  behavior  between  these  two  Mach  number 
regimes. 

The  data  shown  plotted  on  Figure  48  have  been  reduced  by  assuming  that 
the  phase  velocity  equals  the  edge  flow  velocity,  as  is  evident  from  eq.  (13). 
This  is  probably  quite  true  for  the  hypersonic  data,  where  the  prominent 
instability  is  supposed  to  be  of  the  second,  two-dimensional  mode,  and  where 
the  critical  layer  occurs  very  near  the  boundary  layer  edge.  For  the  0  <  M^  < 
3  data,  however,  the  instability  mode  generating  the  data  shown  in  the  figure 
is  thought  to  be  the  first  3-D  mode,  with  phase  velocities  fundamentally 
different  from  the  edge  velocities. 

6.3.3.  The  Amplification  Rates  and  the  Stability  Diagram 
The  amplification  rates 
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were  determined  by  first  curve-fitting  the  A(f;R)  data  vs.  R  by  the  least- 
squares  polynomial  method.  This  curve-fitting  became  a  major  issue  because  of 
our  desire  to  minimize  bias  by  using  the  same  polynomial  degree  for  all  data, 
because  of  the  vide  and  ondulatory  variation  of  the  observed  amplitudes  and 
also  because  of  the  sensitivity  of  the  curve-fits  to  the  number  of  data  points 
included  in  the  fitted  range.  One  can  produce  serious  differences  in  the 
magnitude  of  the  derivatives  of  the  polynomials  A(R;f)  which  enter  (14),  by 
choosing  different  ranges  of  data  to  fit  or  by  choosing  inappropriate 
polynomial  degrees.  Details  of  these  difficulties,  as  they  arose,  are 
described  in  Appendix  D,  and  can  be  appreciated  best  by  scrutinizing  the 
figures  of  that  Appendix  (Figures  D.l,  D.2,  etc.).  The  end  result  was  to 
reduce  data  always  in  the  range  0  <  x  <  10  cm  and  to  use  a  7th  degree 
polynomial  in  that  interval. 

Figures  49  through  51  present  typical  results  comparing  data  points  on 
the  variation  of  A(R;f)  vs.  R  (in  this  case,  A(x;f)  vs.  x)  with  the  polynomial 
chosen  to  represent  each  case.  With  few  exceptions  the  fidelity  of 
representation  is  very  good.  Here  care  must  be  exercised  not  to  compare 
amplitudes  at  different  p^  or  wall  surface  conditions  since  two  different  hot¬ 
wire  probes  were  used.  Although  the  solid  curves  cannot  be  taken  seriously  at 
the  very  extremes  of  the  range  (e.g.  at  x  <  1  cm,  x  >  9  cm)  it  is  already 
clear  from  these  figures  that  the  disturbance  magnitudes  tend  to  increase 
greatly  in  the  downstream  direction  beginning  at  an  x  which  appears  dependent 
on  Pq,  even  though  the  range  0  <  x  <  10  cm  is  almost  always  in  the  laminar 
regime.  Figure  52  is  a  copy  of  the  CRT  screen  display  made  during  the  curve¬ 
fitting  process  and  is  included  to  illustrate  the  ease  with  which  the  operator 


could  examine  the  validity  of  the  curve-fits. 

The  amplification-rate  findings  are  exposed  here  in  a  number  of  ways: 


the  "relief"  plots  of  Figures  53  through  55,  frequency-dependence  plots  of 
Figures  56  and  57,  Reynolds-dependence  plots  of  Figures  58  and  59,  and  finally 
Figure  60  which  attempts  to  extract  the  smooth  vs.  rough  effect. 

To  those  familiar  with  stability  theory,  the  relief  plots  of  Figures  53 
through  55  afford  a  convenient  means  of  judging  the  amplification  rate 
results.  There  are,  in  essence,  three  main  features  worthy  of  note:  first, 
an  up-down-running  ridge  which  is  sometimes  nearly  vertical  (e.g.  p^  *  475 
torr,  smooth  wall)  or  inclined  (p^  •  600  torr,  rough  wall)  in  a  way 
reminiscent  of  the  classic  first-mode  unstable  region  at  low  speeds.  Second, 
there  is  a  weak  evidence  of  a  small  amplification  region  at  low  F  and  low  R, 
which  occasionally  blends  with  the  "principal"  unstable  region  just  mentioned. 
Third,  there  is  a  consistent  and  prominent  onset  of  amplification  at  "large"  R 
which  actually  dominates  the  picture.  Mote  that  this  latter  feature,  marked 
by  a  minimum  (trough)  in  the  amplification  rates,  always  lies  considerably 
upstream  of  the  first  departure  point,  marked  on  Figures  53  -  55  by  a  vertical 
dashed  line  (the  abscissa  on  the  figures  is  obscuring  the  actual 

wetted-length  of  plate  dominated  by  this  feature,  which  length  is  typically 
one-half  of  the  distance  from  the  L.E.  to  the  first  departure). 

Considerable  scrutiny  of  Figures  53-60  is  required  before  the  principal 
lessons  taught  by  them  are  stated;  first,  however,  we  present  on  Figures  60 
and  61  the  stability  diagrams  for  the  smooth  and  rough  walls  in  the  familiar 
R-F  coordinates.  As  the  symbol  key  explains,  these  figures  show  the  points  on 
the  F-R  plane  where  ^  (the  neutral  branches)  with  a  distinction  of  the 

"lower"  and  "upper"  branch  (- increasing  or  decreasing  through  zero, 
respectively).  Dots  and  crosses  identify  points  of  minima  and  maxima  in  -  OCj^, 
respectively;  note  that  a  maximum  in  -oi^may  not  necessarily  lie  in  an 
amplified  region,  and  a  minimum  may  not  necessarily  lie  in  a  damped  region. 
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That  is,  one  can  have  "hilltops  submerKed  under  the  sea"  and  "canyon  bottoms 

on  dry  land."  The  value  of  this  presentation  is  that  the  imminence  of  a 

stable  region  in  the  F-R  plane  can  be  signalled  by  a  minimum  in  -  otXthe  dots) 

long  before  a  neutral  branch  appears.  For  example  at  p^  ~  350  for  the  smooth- 

wall  (Figure  61)  the  lower  neutral  branch  at  R  c::  230  is  signaled  first  by  the 

trough  in  ~  oc  occurring  near  R  *  150,  and  the  second  amplified  region  starting 
i 

at  Rs450  is  signaled  first  by  the  trough  around  R  •=  400. 

It  is  crucial  to  understand  that  the  R  range  shown  on  Figures  61  and  62 
covers  the  region  before  the  first  departure,  i.e.  it  represents  entirely 
laminar,  self-similar  flow.  The  apparent  clutter  of  points  includes  scatter, 
as  usual,  and  probably  some  Instrumentation  problems  at  large  F,  at  the  very 
tops  of  the  figures. 

At  p^  «  600  torr  it  is  clear  that  the  data  clutter  becomes  so  pronounced 
that  it  is  impossible  to  detect  amplified  on  damped  regions  with  any 
certainty.  This  phenomenon  occurs  for  both  the  smooth  and  the  rough  wall.  It 
will  be  recalled  that  at  this  pressure  the  boundary  layer  on  the  sidewalls  was 
completely  turbulent.  The  problem  with  the  600  torr  data  has  therefore  been 
tentatively  ascribed  to  the  increased  disturbance  level  in  the  stream  and,  in 
the  main,  these  data  will  be  excluded  from  further  discussion. 

6.3.4.  The  First  Unstable  Mode  for  the  Smooth  Wall 

We  will  now  attempt  to  identify  and  discuss  specific  features  of  the 
amplification  rate  and  stability  diagram  results  shown  on  Figures  53-62, 
limiting  our  comments  to  the  smooth  wall  unless  otherwise  specified. 

Figure  61  shows  an  unstable  region  in  the  lower  center  of  the  stability 
diagram.  The  neutral  branches,  replotted  on  Figure  63,  are  seen  to  be 
independent  of  unit  Reynolds  number.  Also  plotted  on  the  latter  figure  are 
the  Mg  >  2.2  data  of  Laufer  and  Vrebalovich.  In  their  paper  (Reference  5)  the 


latter  had  compared  their  results  with  early  numerical  results  by  Mack 
(Reference  13)  and  found  good  qualitative  agreement  regarding  the  general 
location  of  this  unstable  region  and  excellent  quantitative  agreement  with 
theory  as  regards  to  the  location  of  the  upper  branch.  Tbeoreticallyi  this 
so-called  first  instability  mode  predicted  by  theory  and  found  experimentally 
by  Laufer  and  Vrebalovich  does  not  shift  greatly  in  the  F-R  plane  as 
changes  from  2.2  to  3;  since  numerical  predictions  for  this  mode  at  ~  3  are 
not  sufficiently  in  existence  from  comparison,  the  similarity  with  the  > 
2.2  data  on  Figure  63  prompt  the  statement  that  the  present  data  plotted  in 
that  figure  represent  the  first  instability  mode. 

An  improvement  over  earlier  measurements  is  noted  in  cni^t  the  minimum 
critical  Reg  (~190)  and  maximum  "amplified  frequency”  F  ( /^  0.000225)  are 
quite  clearly  defined.  Two  fundamental  features  of  the  present  data  deserve 
attention,  however.  One  is  the  shortness  of  the  upper  branch  which  appears  to 
stop  at  Re^  ■  300,  F  ■  0.00013,  thereafter  "merging"  with  the  lower  neutral 
branch  of  a  second  unstable  region  (consult  Figure  61)  to  be  discussed 
shortly.  Another  is  the  presence  of  amplification  at  low  frequencies;  below 
F  *  0.0001  the  present  data  show  amplification  at  all  F  beyond  about  )^80 

whereas  theory  and  Laufer's  tests  indicate  stability,  i.e.  damping,  at  very 
low  F.  As  already  mentioned,  this  low-frequency  ampl if icat ion  bad  already 
been  encountered  at  *  3  by  Kendall,  reportedly  in  an  extreme  form  of 
disturbance  growth  at  all  Reg  .  It  is  not  yet  clear  whether  this  phenomenon 
is  an  aberration  common  to  wind-tunnel  flows  or  a  true  feature  of  stability. 


6.3.5.  The  First-Mode  Amplification  Rates 


Figures  56  and  58  show  typical  amplification  rates  and  Figures  64-65  show 
typical  amplification  spectra;  use  of  R  in  the  abscissa  removes  the  need  to 


choose  between  nominal  and  actual  RCg  .  The  striking  feature  of  Figures  56-58 
is  the  large  rates  observed  near  the  leading  edge  when  “  600  torr,  which 
bears  on  the  discussion  just  concluded  regarding  low-frequency  amplification. 
The  examples  in  the  figures  show  that  at  p^  ”  600  torr, -ol.  is  very  large  near 
the  plate  leading  edge.  This  would  indicate  an  increased  efficiency  of  the 
leading  edge  region  as  a  conduit  and  cultivator  of  spurious  noise  into  the 
boundary  layer,  when  the  latter  is  irradiated  by  external  disturbances. 

An  attempt  has  been  made  in  Figure  66  to  compare  the  measured 
amplification  rates  with  theory.  The  only  numerical  amplification  rate 
results  available  for  the  first,  three-dimensional  mode  at  Mach  3  are  due  to 
L.  Mack,*  and  they  are  shown  on  Figure  66  also.  The  computation  pertains  to 
oblique  waves  in  the  range  55*  <  65°,  for  which  the  rates  are  presumably 

the  highest.  At  the  lower  R  (e.g.  R  ■  300)  the  present  data  are  obviously 
influenced  by  the  low-frequency  gain  phenomenon,  while  at  higher  R  (e.g. 
R  ■=  500)  there  is  an  apparent  influence  of  other  amplification  mechanisms,  to 
be  discussed  later.  Otherwise,  however,  the  data  of  Figure  66  are  not  far 
from  the  theory,  thus  justifying  their  identification  with  the  first  mode. 

The  practice  of  plotting  neutral  branches  alone  does  not  yield  the 
maximum  information  on  stable  and  unstable  regions.  Additional  knowledge  can 
be  gained  from  other  topographical  features  of  the  stability  diagram  by 
including  the  extrema  in  the  amplification  rates  (the  "ridge  tops"  and  "canyon 
bottoms").  A  case  in  point  is  Figure  67  where  the  maxima  in  — o<.(the  maxima 
found  in  Figure  61  for  example)  are  shown.  The  fact  that  such  maxima  usually 
run  down  the  middle  of  unstable  regions  is  confirmed  by  the  cluster  of  points 
"down  the  middle"  of  the  unstable  region.  The  Figure  also  shows  the 
topographical  relation  among  this  maximum  amplification  line,  the  maximum 

*  The  author  is  indebted  to  Mr.  Mack  for  these  unpublished  computations. 


amplitude  line  (see  Figure  47)  and  the  neutral  branches;  it  also  finds  little 
differences  between  use  of  the  actual  and  nominal  Re0  . 

The  important  feature  of  Figure  67  is  that  the  p^  “  600  torr  data  (Re'  * 
56100/cm),  despite  scatter,  are  also  well  represented  along  the  line  of 
maximum  -  oc^.  Thus,  whereas  neutral  branches  at  this  Re'  are  hard  to  discern 
(see  Figure  61),  a  strong  hint  toward  the  presence  of  the  first  mode  is 
afforded  by  this  type  of  plot  even  in  this  case  of  the  boundary  layer 
irradiated  by  stream  noise. 

6 .3 .6 .  Effect  of  Roughness .  Turbulence  and  Unit  Rovnolds  No.  on  the  First  Mode 

Reviewing  the  effect  of  unit  Reynolds  number  Re'  on  the  first  mode,  we 
realize  that  we  can  at  best  study  this  effect  by  comparing  the  Re'  29400/ cm 
and  43900/cm  data;  use  of  the  56100/cm  data  is  questionable  because  the  latter 
also  includes  the  stream-irradiation  effect  (the  "turbulence").  Within  this 
context,  we  have  seen  that  Re'  changes  do  not  affect  the  location  of  the 
unstable  region  on  the  stability  diagram.  There  appears  to  be,  however,  an 
effect  on  the  amplification  rates.  From  Figures  64  and  65  we  see  that  the 
higher  Re'  data  are  more  damped  especially  at  low  R  or  Re..  This  must  be 
connected  with  the  observed  effect  of  Re'  on  transition  (Figure  16). 

For  this  particular  experiment  it  appears  that  noise  irradiation 
increases  the  "noise"  of  the  measurement,  causes  amplification  near  the 
leading  edge  at  practically  all  frequencies  and  obscures  the  neutral  branch 
location.  The  end  effect,  however,  is  not  to  aid  and  abet  the  first  mode  in 
accelerating  the  transition  process;  the  transition  data  of  Figure  16  do  not 
show  any  abrupt  changes  of  behavior  when  sidewall  radiation  sets  in. 
Furthermore,  the  stream  disturbance  levels  of  Figure  32  show  only  a  moderate 
increase  in  level  when  sidewall  radiation  arrives  at  the  point  of  measurement. 


By  contrast,  the  data  of  Beckwith  and  Kendall  (Reference  24)  show  a  much 
larger  stream  noise  increase  upon  irradiation,  accompanied  by  markedly  earlier 
transition.  At  this  point,  therefore,  the  stream  noise  effect  in  the  present 
test  appears  much  more  benign  than  expected. 

The  rough  wall  results  are  much  less  conclusive  than  Figure  62  indicates. 
Since  the  roughened  surface  did  not  start  until  3  cm  from  the  leading  edge  (R 
^300  at  Pq  =  350  torr),  one  is  at  a  loss  to  explain  the  unstable  region  shown 
on  Figure  62.  The  same  puzzlement  can  be  expressed  regarding  the  elevated 
amplification  rates  over  the  rough  wall  at  low  R,  as  per  the  amplification 
spectra  of  Figures  64  and  65.  Discussion  of  the  effect  of  roughness  on 
instability  will  therefore  be  postponed  until  further  data  are  discussed 
later. 


6.3.7.  The  Second  Instability  Region 

It  is  evident  from  the  overall  stability  diagrams  of  Figures  61  and  62 
that  a  second  amplification  region  exists  to  the  right  of  the  first  mode 
described  in  the  previous  section.  We  will  refer  to  this  additional  region  of 
amplification  as  the  "second  instability  mode"  partly  for  distinguishing  it 
from  the  first  mode  results,  and  partly  because  of  evidence  of  its  connection 
to  the  second  mode  of  instability  in  hydrodynamic  theory. 

Evidence  of  the  second  mode,  as  R  increases,  first  appears  as  the  locus 
of  minima  in -oc.located  nearly  vertically  in  Figures  61  and  62  at  R~400, 
followed  by  a  lower  neutral  branch  in  the  region  450  <  R  <  550.  The  latter 
branch  has  been  for  clarity  replotted  on  Figure  68  (and  also  on  Figure  69  for 
the  rough  wall).  It  is  seen  that  this  amplification  phenomenon  considerably 
precedes  the  first  departure,  that  it  "intercepts"  the  upper  neutral  branch  of 
the  first  mode,  that  it  extends  to  high  frequencies  (of  order  F  ^  0.0004)  and 
that  it  is  strongly  Reynolds-number  dependent. 


Typical  amplitude  and  amplification  data,  which  generated  the  stability- 
diagram  plots  of  Figure  61,  are  shown  on  Figures  70  through  73.  The 
amplitudes  A(f)  shown  on  Figure  70  for  typical  cases  at  p^  *  350  and  475  torr 
are  especially  illuminating  because  they  show  the  relative  insignificance  of 
the  first  compared  to  the  second  mode.  The  former  creates  a  rather  minor 
disturbance  in  the  evolution  of  fluctuations  along  the  plate,  while  the  latter 
involves  a  pronounced  increase  in  the  signal.  Note  that  the  wetted  length 
of  plate  over  which  the  second  mode  is  active  in  the  self-similar  laminar 
region  is  quite  long.  For  example.  Figure  70  shows  that  for  F  =  0.000207  the 
first  mode  activity  lies  between  R  =  260  and  330,  roughly,  corresponding 
to  1.4  cm  at  p^  =  350  torr;  while  the  second  mode  is  active  from  R  “  460 

and  is  still  so  at  the  first  departure  (R  ^  600),  for  an  "active  length"  ZkX  = 
5  cm.  When  p^  =  475  torr,  from  Figure  70  we  get  Ax  =  2.6  cm  for  the  region 
of  first  mode  activity  in  300  <  R  <  450,  and  Ax*  3  cm  for  the  second  mode 
activity  region  in  525  <  R  <  640.  Similar  plots,  stressing  the  relative 
extent  of  the  two  modes  in  the  physical  coordinate  x  along  the  plate  are  shown 
on  Figures  71-73.  The  top  plot  of  each  figure  is  one  of  amplitude  and  the 
bottom  two  of  amplification.  The  two  "bulges"  in  the  latter  are  the  two 
instability  modes  under  discussion.  The  predominance  of  the  second  mode  over 
the  first  is  evident  even  for  F  as  low  as  0.000056.  Note,  too,  that  the 
second  mode  becomes  more  prominent  as  F  increases.  For  example  at  F  *  0.00021 
which  is  near  the  highest  amplified  frequency  for  the  first  mode,  considerable 
amplification  of  the  natural  disturbances  due  to  the  second  mode  alone  occurs 
in  the  range  7  <  x  <  12  cm  when  Re'  “  29400/cm.  It  is  also  evident  that  at 
very  low  frequencies,  say  F  <  0.00005,  the  two  modes  merge  and  amplification 
(starting  at  an  unexpectedly  early  point,  as  already  discussed)  continues 
uninterruptedly  throughout  the  self-similar  region  since  the  second  mode  "sets 
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in"  before  the  first  mode  is  done  with  amplification.  This  merging  process 
makes  it  difficult  to  discuss  any  further  the  topography  of  the  second  mode  at 
the  low  frequencies,  a  topic  which  will  be  left  aside  until  total 
amplification  is  discussed  later. 

According  to  the  figures  shown  so  far,  the  amplification  rates  for  the 
second  mode  appear  to  be  much  higher  than  those  of  the  first  mode,  especially 
at  large  F  and  R.  To  show  numerical  results,  the  amplification  spectra  of 
Figures  64,  65,  etc.  have  been  extended  by  plotting  on  Figures  74  to  77  such 
spectra  to  R  =  600.  For  these  plots,  the  computations  were  made  using  Option 
11  in  the  program  STABLE02  (see  Appendix  C)  but,  in  contrast  to  the  first  mode 
presentations  of  Figures  64,  65,  etc.,  the  curve-fitted  data  extended  to  x  = 
15  cm  (for  Figures  64,  65,  etc.  the  range  extended  only  to  10  cm,  which  is 
where  the  second  mode  becomes  most  active  and  might  have,  therefore,  missed  or 
misrepresented  the  second-mode  results). 

Figures  74  through  77  show  the  second-mode  amplification  spectra  —  vs. 
F  in  plots  similar  to  those  of  Figures  64,  65,  etc.  The  frequency  range 
extends  to  F  =  Q0004  and  the  R  values  displayed  are  R  =  300,  400,  500  and  600. 
These  ranges  again  display  some  of  the  amplification-spectrum  results  shown 
earlier  but  in  a  way  which  now  allows  a  good  comparison  between  the  two  modes, 
and  also  of  the  individual  mode  behavior  along  the  plate.  If,  for  the  moment, 
we  concentrate  on  the  smooth  wall  results  (Figures  74  and  75)  we  note  the  ease 
with  which  the  first  mode  shift  toward  the  lower  frequencies  can  be  followed 
as  R  increases.  Maximum  rates  for  this  mode  reach  about  0.005.  The  second 
mode  appears  to  set  in  almost  abruptly  at  R  —  450  and  to  involve  maximum  rates 
of  order  0.012.  Therefore,  the  second  mode  predominates  at  locations  on  the 
model  ^ust  preceding  the  first  departure  (transition),  and  involves  much 
higher  amplifications  than  ever  attained  by  the  first  mode.  On  occasion  (e.g. 
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Figure  75  at  R  =  600  for  Re'  =  43900/cni,  which  occurs  just  before  the  first 
departure)  the  first  mode  even  appears  as  a  minor  superposition  on  the  second 
mode  spectrum.  This  is  very  significant  because  it  implies  that  the  second 
mode  activity  extends  to  low  frequencies  (lower  than  F  =  0.0001)  where,  as  we 
shall  later  see,  one  must  look  for  the  causes  of  transition.  In  other  words, 
the  second  mode  may  here  be  at  the  least  a  participant  in,  if  not  the  cause 
of,  transition. 

A  final  point  concerns  the  chosen  "first  mode"  and  "second  mode" 


terminology.  In  the  context  of  stability  theory  modes  are  specific  "groups" 
of  eigenvalue  solutions  of  the  small  disturbance  equations  which  can  coexist 
at  a  certain  R  with  different  amplification  behavior;  furthermore,  at  the  same 
R  the  amplification  spectrum  (  —  ,  F)  of  a  mode  can  depend  on  the 

disturbance  wavefront  orientation  relative  to  the  flow  vector.  It  is 
experimentally  difficult  to  attribute  an  observed  disturbance  amplification  to 
the  action  of  any  particular  mode,  although  sometimes  the  observed 
amplification  spectra  can  be  decomposed  into  individual  mode  contributions.  A 
familiar  case  in  point  is  the  plot  on  Figure  77,  taken  from  Reference  9, 
showing  amplification  spectra  for  =  6.8  to  7.7  and  R  -  1740.  For  reasons 
clear  from  this  plot,  the  community  has  interpreted  the  data  shown  on  it  as  a 
superposition  of  two  theoretical  instability  modes,  a  three-dimensional  one 
with  orientation  =  55  and  a  two-dimensional  one  with  'p  ’=  0  .  Based  on 
obvious  similarities  between  this  plot  and  the  present  amplification  spectra 
from  Figures  74  through  77,  it  is  evident  that  the  latter  also  shows  the 
merging  of  two  different  modes,  shown  by  single  or  double  arrows  and  called 
"first"  and  "second"  respectively.  From  what  has  been  shown,  it  is  also 
evident  that  the  use  of  the  word  "first"  (for  "first  3-D")  is  justified,  but 
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we  know  little  at  this  point  on  whether  the  second  mode  is  two-  or  three- 
dimensional.* 

6.3.8.  The  Character  and  Mach-Number  Dependence  of  the  Instability 

It  is  seen  from  the  previous  discussion  that  the  unstable  region  found  is 
rather  complex,  consisting  of  two  merged  instability  modes  which  provide  two 
separate  peaks  in  the  amplification  spectra.  At  lower  R  (say  R  <  400)  these 
two  modes  are  separated  by  a  damped  region  which  disappears  when  R  increases. 

As  was  already  indicated  by  the  example  of  Figure  77,  this  type  of  mixed¬ 
mode  or  merged-mode  instability  is  common  at  hypersonic  Mach  numbers.  The 
experiments  of  this  author  (References  6,  7  and  8)  and  of  Stetson  et  al. 
(Reference  9)  show  the  existence  of  a  large  region  of  instability  in  the  range 
1000  <  R  <  2000  in  which  the  first  and  second  modes  co-exist.**  The  present 
day  experimental  knowledge  of  the  Mach  number  dependence  of  tbe  stability 
diagram,  typified  by  the  present  data  and  those  of  References  3,  5,  8  and  9  is 
pictured  on  Figures  78  and  79.  The  interesting  feature  here  is  that  the 
"second  mode"  lower  neutral  branch  claimed  by  the  present  data  appears  to  form 
the  natural  extension  of  the  hypersonic  results  (at  present  limited  to  R  > 
1000)  at  the  lower  Reynolds  numbers. 

*There  are,  strictly  speaking,  several  possibilities  for  second  mode 
mechanisms  including  some  not  related  to  Tollmien-Schlichting  phenomena  such 
as  three-dimensional  or  Tay lor-Goertler  instabilities.  There  may  also  be 
tunnel-peculiar  phenomena,  such  as  finite-span  problems  wherein  instabilities 
generated  by  three-dimensionalities  at  the  plate-sidewall  intersection 
propagate  to  the  plate  centerline. 

**The  authors  of  these  references  found  additional  unstable  regions,  consisting 
of  "higher  harmonics"  of  the  main  instability  under  discussion  here. 


This  opinion  is  reinforced  by  Figure  79  where  the  abscissa  is  the 
Reynolds  number  based  on  the  kinematic  thickness  of  the  boundary  layer.  This 
plot  was  motivated  by  Laufer's  original  observation  that  such  scaling  brings 
the  neutral-branch  data  of  References  3  and  5  into  coincidence.  To  extend 
this  test  to  the  present  data,  the  latter  were  converted  from  R-values  to  Re^ 
values  using  the  measured  thickness  S  .  The  hypersonic  results  were 
converted  to  the  coordinate  Re  by  using  this  author's  measured  relation  (see 
Reference  8,  Figure  3)  between  Reg  and  R; 

Reg  sil2R 

Figure  79  thus  shows  the  neutral  branch  location,  on  the  F-Reg  plane,  of  the 
instabilities  most  likely  to  be  experimentally  detected  over  a  wide  range  of 
Mg.  Of  course,  observability  is  mainly  a  question  of  amplification  rate;  low- 
speed  measurements  cannot  discern  the  '*hyper sonic”  instability  because  the 
latter  has  no  amplification  (is  not  present)  at  low  Mg. 

6.3.9.  Roughness  and  Reynolds  Number  Effects  on  the  Second  Mode 


It  has  already  been  mentioned  that  the  lower  neutral  branch  of  the  second 
mode  is  Reynolds-number  sensitive.  The  results  of  Figure  68  supporting  this 
claim  pertain  to  the  smooth  wall,  and  in  this  case  the  effect  is  a  systematic 
shift  of  the  neutral  curve  toward  larger  R  as  Re'  increases.  It  is  worth 
repeating  here  that  this  constitutes  a  possible  explanation  of  the  unit 
Reynolds  number  effect  on  transition. 

The  corresponding  effect  of  Re'  for  the  rough  wall  gives  the  surprising 
results  of  Figure  69.  Increasing  Re'  in  this  case  changes  completely  the 
character  of  the  lower  neutral  branch,  from  one  with  positive  slope  in  the  F, 
R  plane  when  Re'  =  29400/cm,  to  one  with  negative  slope  when  Re'  ■=  43900/cm. 
No  explanation  is  available  at  this  stage,  and  no  definite  conclusions  were 
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therefore  pursued  about  the  effect  of  Re'  when  the  wall  was  rough. 

6 .4 .  Boundarv-Laver  Response  and  the  Transition  Question 

While  the  presentation  and  discussion  of  the  data  have  so  far  dealt  with 
details  of  the  small  disturbance  development  along  the  laminar  boundary  layer, 
we  are  ultimately  interested  in  the  cumulative  effect  of  amplification  or 
damping  and  in  the  total  amplification  experienced  by  each  wave  number  between 
entry  into  the  boundary  layer  and  arrival  at  the  first  departure  zone.  This 
net  cumulative  effect,  presently  termed  the  "boundary  layer  response", 
associates  the  total  "gain"  experienced  by  each  wave  number  with  its 
likelihood  of  triggering  transition,  the  presumption  being  that  those 
disturbances  displaying  the  maximum  total  gain  before  the  first  departure  are 
responsible  for  causing  the  transition  phenomenon.  Typical  targets  for  such 
gain  computations  would  be  to  find  a  net  factor  by  which  a  disturbance  must 
increase  before  causing  transition  or  even  an  absolute  disturbance  magnitude 
needed  to  cause  transition. 

The  disturbance  amplitude  gain  within  the  boundary  layer  can  be  discussed 
with  the  aid  of  Figure  80.  The  line  of  constant  y/S  along  which  data  were 
obtained  can  be  roughly  approximated  with  a  "disturbance  track"  originating  at 
a  point  X  =  Xjj  (entrainment  point)  on  the  layer  edge,  close  to  the  plate 
leading  edge,  where  imall  convected  or  radiated  stream  disturbances  enter  the 
boundary  layer.  Comparison  of  the  disturbance  amplitude  at  any  point  on  this 
track  with  its  "entry"  amplitude  at  x^  will  give  the  cumulative  effect  of 
boundary  layer  instability  on  the  disturbance  up  to  that  point;  thus  the  gain 
up  to  a  certain  point  for  a  Fourier  component  F  is  simply  the 

"spectrum  ratio"  A(R;F)  to  A(Rq;F): 


G(R;F) 


A(R;F) 

A(Ro;F) 


;  Ro3R(x*=Xq) 


(16) 


Since  the  entrainment  point  is  difficult  to  fix  or  define,  gain 
calculations  for  the  present  data  were  made  by  picking,  for  x^,  a  point  x^ 
near  the  most  forward  probe  position;  for  convenience,  the  x^  was  chosen  so 
that  Rq  =  150  always  (this  gave  x^  “  0.8  cm  for  Re'  =  29400/cm  and  0.8  cm  for 
Re'  =  43900/cm).  Results  are  shown  on  Figures  81  and  82.  As  has  been 
increasingly  evident  in  this  report,  the  disturbances  have  gained  little  in 
amplitude  between  the  leading  edge  and  R  =  450.  For  the  smooth  wall  at  Re'  * 
43900/ cm  all  disturbances  have  in  fact  experienced  net  damping,  up  until  R  “ 
450  (Figure  81).  The  disturbance  amplitude  gain  seems  to  be  large  only  when 
the  second  mode  becomes  active,  in  the  interval  between  R  -450  and  the  first 
departure.  Figures  81  and  82  show  that  the  maximum  gain  at  the  first 
departure  is  of  order  2.5  to  3  and  is  confined  to  the  low  frequencies  (F'^ 
0.000G5).  This  result  is  of  course  obvious  even  at  the  early  stages  of  this 
narrative  during  discussion  of  the  amplitude  histories  (Section  6.3.3  and 
Figures  49,  50,  etc.). 

The  factor  of  3  is  too  small  to  be  a  credible  mechanism  for  initiating 
transition  in  the  boundary  layer.  It  must  be  observed  from  the  data  shown  so 
far  that  the  disturbance  magnitude  at  the  beginning  of  the  process  (i.e.  near 
Xq  or  x^,  that  is  near  the  leading  edge)  is  already  substantial  compared  to 
the  stream  disturbance  level.  Previous  observations  (Reference  12)  have 
already  shown  that  much  larger  gain  factors  obtain  when  the  disturbance 
magnitudes  just  before  transition  are  compared  with  the  stream  disturbance 
level  Ag,  which  must  be  correctly  considered  to  be  the  "initial"  disturbance 
magnitude.  Mack  (Reference  14)  has  formulated  a  "forcing"  variant  of 
stability  theory  which  accounts  for  the  boundary  layer  stability  together  with 
the  stream  noise  environment.  At  a  Mach  number  of  4.5,  Kendall  has 
qualitatively  verified  the  forcing  theory  and  found  gains  A/Ag  just  before 


transition  which  amount  to  factors  10-100  depending  on  frequency. 

To  investigate  the  total  gain  in  the  boundary  layer  of  a  disturbance 
which  starts  out  with  a  magnitude  A^Cf)  in  the  stream,  we  present  the  spectra 
A(F)/Ag(F)  of  Figures  83  and  84.  The  gain  intrinsic  to  the  boundary  layer, 
just  discussed  above,  can  be  seen  in  these  Figures  as  the  difference  between 
A/Ag  at,  say,  R  200  and  the  first  departure  R  (R  ^  600  for  smooth  wall). 
For  example,  for  the  smooth  wall  at  Re'  “  29400/cm  the  gain  at  R  *  229,  F  * 
0.00005  is  about  5  and  at  R  ”  592  it  is  15,  giving  the  factor  3  previously 
discussed.  The  principal  message  of  Figure  83  is  that  at  low  R,  of  order  200 
(or  of  order  150  in  the  vicinity  of  x^.  Figure  84)  the  disturbance  magnitude 
is  a  factor  of  as  high  as  5  larger  than  that  in  the  stream.  Note  that  the 
growth  is  nil  or  "sluggish"  for  a  certain  range  of  R  thereafter.  Translating 
R  into  X  for  the  typical  case  of  Re'  =  29400/ cm,  we  can  put  together  the 
following  physical  picture  of  disturbance  gain  for  low  F:  upon  "entering"  the 
layer  between  the  L.E.  and  x  •=  1.8  cm  (R  229)  the  disturbance  increases  by  a 
factor  of  5  in  size,  and  gains  an  additional  factor  of  3  or  so  between  x  *  1.8 
cm  and  x  *  12  cm  (R  592). 

This  point  is  further  demonstrated  in  the  cross-plot  A/A^  vs.  R  of  Figure 
84  for  F  “  0.00005.  Qualitatively  similar  results,  plotted  on  the  same 
Figure,  have  been  reported  by  Kendall  (Reference  12).  The  latter's  data  show 
a  slower  gain  process,  with  "non-similarity"  at  Mg  =  3  and  Re'  *  43,000/cm 
indicated  at  R  900  vs.  our  own  first  departure  point  around  R  •*  600.  For 
this  reason  Kendall's  final  gain  before  transition  seems  to  be  A/Ag  25  at 
this  F,  while  ours  is  about  15.  In  other  words,  Kendall's  finding  says  that 
transition  begins  around  19  cm,  after  the  disturbance  magnitude  has  increased 
25-fold,  while  we  find  it  occurring  at  12  cm  following  a  15-fold  increase. 
This  difference  might  be  attributed  to  differences  between  the  two  wind- 
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tunnels,  although  in  both  cases  there  vas  no  irradiation  from  turbulent 
sidewall  boundary  layers.  On  the  other  hand,  there  may  be  other,  simpler 
explanations  dealing  with  interpretation  of  the  measurements.  For  instance, 
it  is  not  clear  from  Kendall's  article  that  his  "non-similarity"  corresponds 
to  our  "first  departure."  More  important,  in  both  cases  the  changes  in  hot¬ 
wire  sensitivity  and  frequency  response  are  not  really  accounted  for  in 
measuring  A  and  A^;  strictly  speaking,  therefore,  all  data  on  Figure  84  are 
qualitative.  Finally,  whereas  our  A^  was  measured  just  over  the  plate 
boundary  layer,  Kendall  himself  cautions  that  his  A^  vas  measured  ahead  of  the 
plate  and  thus  is  probably  overestimated;  in  comparing  his  data  with  theory  he 
increases  his  A/A^  by  a  factor  of  2  for  comparison  purposes  (Reference  12,  p. 
245).  If  the  solid  curve  is  also  increased  by  2  on  Figure  84,  then  Kendall's 
and  the  present  data  come  into  much  better  agreement. 

No  theoretical  results  from  Mack's  combined  forcing-and-stability  theory 
are  available  for  M^  =  3^  but  a  comparison  of  the  present  data  with  his  theory 
(quoted  by  Kendall,  Reference  12,  p.  295)  at  M^  **  4.5  is  also  shown  on 
Figure  84,  as  a  matter  of  some  interest. 

Conclusions 

Based  on  the  results  obtained  from  the  present  measurements,  the 
following  conclusions  can  be  drawn  for  the  flat-plate  boundary  layer  at 
Mach  3: 

1.  An  unstable  region  has  been  found  which  has  been  identified  with  the  first 
unstable  mode.  Supporting  evidence  consists  of  the  shape  and  location  of 
this  region  on  the  stability  diagram  and  the  agreement  of  its 
amplification  rates  with  the  available  theoretical  expectations.  Minimum 
critical  Re.  for  this  region  is  190  and  maximum  amplified  frequency  is  F  ■= 
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0.000225,  independently  of  unit  Reynolds  number. 

2.  An  extensive  second  unstable  region  was  also  found  with  amplification 
rates  much  larger  than  those  of  the  first  mode,  extending  to  frequencies 
higher  than  F  =  0.00035  and  strongly  dependent  on  unit  Reynolds  number. 
This  region  separates  from  the  first-mode  at  low  R,  merges  with  the  latter 
at  high  R  and  appears  to  connect  with,  and  forms  the  low-R  terminus  of, 
the  merged  first-  and  second-mode  instability  found  earlier  at  high  Mach 
numbers.  For  this  reason,  this  second  unstable  region  has  been  associated 
with  the  second  instability  mode. 

3.  The  net  gain  in  disturbance  amplitude  from  the  free  stream  to  the  first 
departure  is  about  15  and  is  confined  to  a  frequency  band  centered  near  F 

0.00005.  Of  this  gain,  a  factor  of  about  3  is  due  to  the  detected 
boundary  layer  instabilities  and  is  almost  entirely  due  to  the  second 
mode;  the  remaining  factor  of  5  is  experienced  very  near  the  leading  edge 
of  the  plate. 

4.  Transition  delay,  expressed  as  an  increased  Re^  of  transition  as  the  unit 
Reynolds  number  increases,  seems  to  be  caused  mostly  by  a  shift  to  higher 
Reynolds  number  of  the  second  unstable  mode  and  much  less  so  by  decreases 
in  the  first-mode  amplification  rates  as  the  unit  Reynolds  number 
increases. 

5.  When  computed  with  the  stream  rather  than  the  phase  velocity,  the  maximum- 
amplitude  waves  ("laminar  waves"  or  '*rollmien-Schlichting  waves")  of  the 
first  mode  had  a  length  of  about  30S  in  agreement  with  the  previous 
stability  data  for  <  3;  by  contrast  the  second-mode  wavelengths, 
reported  so  far  only  at  hypersonic  speeds,  have  a  primary  wavelength  of 
about  25. 

6.  The  effect  of  sound  irradiation  from  turbulent  sidewall  boundary  layers 
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was  a  doubling  of  the  stream  fluctuations  without  any  discernible  effect 
on  transition  on  the  model.  In  the  boundary  layer  sucb  irradiation  caused 
large  increases  in  the  amplification  rates  near  the  leading  edge,  and 
scatter  of  data  which  tended  to  obscure  the  first  mode,  but  left  the 
second  mode  largely  unaffected. 

A  clearly  recognizeable  effect  of  the  surface  roughness  was  to  make  the 
layer  more  selective,  i.e.  to  produce  peaks  in  the  fluctuation  spectra. 
Otherwise,  peculiar  behavior  observed  of  the  boundary  layer  stability  with 
roughness,  and  its  connection  to  the  observed  earlier  transition  has  been 
left  uninterpreted  because  of  the  unique  (non-general izable)  nature  of  the 
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Figure  2.  Axisymmetric  model  installed 
in  the  SWT  test  section. 
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Figure  3.  Axisymmetric  model,  shown 
disassembled,  with  a  variety  of  rough 
surface  afterbodies. 
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Figure  14,  Typical  velocity  profiles 
(left)  and  Reg  variation  (right)  for  the 
rough  wall. 
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Figure  17.  Transition  data  on  the  flat 
plate  superposed  on  a  correlation  by 
Pate  (Reference  23,  p.  258)  for  various 
wind-tunnels. 
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Figure  20.  Unit  Reynolds  number  effect 
on  the  flat  plate  transition  superposed 
on  a  correlation  by  Beckwith  et  al. 
(Reference  24,  Figure  17), 
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Figure  22.  Basic  elements  in  the  pitot- 
probe  location  on  the  velocity  profile. 
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Figure  25,  Effect  of  roughness  on  the 
measured  integral  properties  at  350  torr 
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Figure  29.  Schematic  of  the  hot-wire 
anemometer  probe. 
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Figure  30.  The  test  set-up  for  the  hot¬ 
wire  measurements 


Unresolved  wideband  hot-wire 
OSS  the  laminar  boundary 
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Figure  32.  Reynolds  number  effect  on 
stream  noise  detected  just  outside  the 
boundary  layer. 
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Figure  34.  Wideband  hot-wire  output 
variation  along  plate  (rough  wall). 
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Figure  35.  Reynolds  number  effects  on 
wideband  signal  along  plate  (smooth  wall). 
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Figure  36.  Reynolds  number  effects  on 
wideband  signal  along  plate  (rough  wall) 
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Figure  38.  As  in  Figure  37;  smooth  wall 
at  475  torr. 
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Figure  43.  Typical  spectra  in  the 
laminar  boundary  layer;  smooth  wall 
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Figure  48.  Prominent  "T-S"  wavelengths  X  visible  in  boundary 

layers  (here  computed  using  the  stream  as  opposed  to  the 
phase  velocity)  group  themselves  in  clusters  depending  on 
the  active  instabil ity  mode.  "Higher  harmonic"  cluster 
found  at  hypersonic  speeds  by  Demetriades  and  later  by 
Stetson  et  al  with  A  <  S  ,  lie  below  the  bottom  margin. 
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Figure  49.  Typical  curve-fits  of  A(x;f) 
at  350  torr. 
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Figure  54.  As  in  Figure  53  for  475  torr 
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Figure  56.  Smooth-wall  amplification- 
rate  dependence  on  R  for  various  frequen 
cies. 
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Figure  57.  Rough-wall  amplification-rate 
dependence  on  R  for  various  frequencies. 
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Figure  62.  Unedited  stability  diagram 
for  the  rough  wall. 
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Figure  66.  First-mode  amplification-spectrum  data  compared  with 
numerical  results  ("theory")  supplied  by  Mack.  Low- 
frequency  anomaly  has  been  also  previously  reported. 
High-R  spectrum  reflects  contamination  by  the  second 
mode.  118 
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Figure  67.  Maximum  amplification-rate 
points,  smooth  wall. 
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stories  illustrate  the  minor  role  played  by 
de.  Note  position  where  boundary  layer  velo' 
'St  departs"  from  the  Blasius  theory. 


Figure  71.  Typical  smooth-wall  amplitude 
(top)  and  amplification-rate  histories 
(middle  and  bottom) ,  350  torr. 
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Figure  73.  As  in  Figure  71,  600  torr 
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Figure  79.  Same  as  Figure  78  but  plotted  vs  Re*  as  suggested  by  Laufer 
to  bring  all  low-M  data  into  coincidence. 


A CR- 150) -AMPLITUDE  AT  R-150 
<X-0.  8  CM) 


_j - 1 - 1 - 1 - 1 — 

-5  1  1.5  2  2.5 

XIO'^^-F 

Figure  82.  Boundary  layer  disturbance 
amplitude  gain  referred  to  the  amplitude 
at  R  =  150,  rough  wall 


Figure  83.  Disturbance  amplitude  gains 
referred  to  the  stream  disturbance  levels 


APPENDIX  A 


Method  of  Decreasing  the  Model-Wall  Interference 

Boundary  layers  tend  to  separate  in  an  adverse  pressure  gradient*  A 
classical  example  of  separation  can  be  found  around  the  base  of  an  object 
protruding  from  an  otherwise  flat  surface  over  which  the  boundary  layer  forms. 
In  the  present  experiment,  the  plate  model  spanned  the  test  section  from  one 
sidewall  to  the  other.  The  mechanism  mentioned  above,  which  is  potentially 
capable  of  separating  the  layer,  is  present  at  the  two  points  where  the  plate 
leading  edge  meets  the  two  sidewalls. 

Observation  made  during  this  experiment  showed  anomalous  flow  along  the 
model  centerline  beginning  at  about  x  6  cm  and  present  at  the  two  lower 
pressure  conditions  (P^  =  475,  350  mm  Hg  abs.).  Mach  wave  tracing  indicated 
that  disturbances  at  that  location  would  originate  at  the  points  of  model-wall 
intersection.  Furthermore,  it  has  been  established  (Reference  15)  that  at  the 
pressures  quoted  above,  the  sidewall  boundary  layer  is  largely  laminar  whereas 
at  P^  =  600  mm  Hg  it  is  turbulent.  This  gave  added  strength  to  the  view  that 
the  anomaly  was  caused  by  sidewall  boundary  layer  separation  at  the  model-wall 
intersection  points  since  turbulent  (rather  than  laminar)  boundary  layers  are 
much  more  resistant  to  separation.  A  schematic  of  this  situation  appears  on 
Figure  A.l. 

According  to  the  above  remarks,  the  diagnosis  of  the  problem  also 
suggests  its  cure  which  would  be  to  cause  the  sidewall  boundary  layer  to 
become  turbulent  upstream  of  the  model-wall  intersection.  This  can  be 
accomplished  by  "tripping"  the  boundary  layer  by  means,  for  example,  of 
artificial  roughness.  Fortunately,  a  great  deal  of  information  about 
transition-tripping  of  the  test-section  boundary  layers  in  the  SWT  was 
available,  and  specifically  about  using  surface  roughness  for  the  purpose 


(Reference  17).  A  0.5"  x  4"  strip  of  120-grit  sandpaper  was  therefore  affixed 
on  the  glass  sidewall  by  adhesive  with  the  leading  edge  of  the  strip  upstream 
of  and  the  trailing  edge  downstream  of  the  sonic  throat;  a  similar  strip  was 
placed  at  the  corresponding  position  on  the  opposite  sidewall.  According  to 
Reference  17,  the  streamwise  extent  of  these  strips  should  be  sufficient  to 
generate  turbulence  downstream  of  them  all  the  way  to  the  diffuser  entrance 
along  a  narrow  band  which  covers  the  model  wall  intersection  point.  The 
placement  at  the  sonic  throat,  where  the  nozzle  boundary  layers  are  thin, 
allows  for  the  use  of  fine-grit  sandpaper  (grit  finer  than  120  could,  in  fact, 
be  used  according  to  Reference  17)  and  thus  minimizes  wavelets  radiated  into 
the  flow  by  the  grit  particles.  Futhermore,  the  leading  edge  of  the  strip, 
being  in  the  subsonic  portion  of  the  nozzle,  generates  no  shock  waves,  and  the 
strips  themselves  were  in  anyway  so  far  upstream  of  the  nozzle  that  they 
caused  no  inviscid  disturbances  to  the  model  flow. 

Subsequent  testing  showed  that  the  installation  of  these  trips  restored 
the  measured  boundary  layer  growth  on  the  model  to  that  expected  from  theory. 
All  data  shown  in  this  report  were  taken  with  the  trips  in  place. 


APPENDIX  B 


Mean  Flovfield  Computer  Programs 

In  order  of  presentation,  the  following  are  included  here: 

B.l  The  PITOXXXX  data  file  (example),  which  is  used  as  an  input  to  the  LAMBL- 
series  programs.  Note:  This  file  is  translated  into  an  indentical  one 
put  without  line  numbers,  by  the  DATAFORM  program,  before  insertion  in 
such  programs. 

B.2  The  LAMBL2  program.  When  asked  to  RUN,  this  program  requests  a  PITOXXXX 
input.  The  program  produces  full  printouts  of  the  boundary  layer 
profile,  using  offset-type  data  reduction;  a  typical  output  is  shown 
under  B.4.  (This  program,  as  well  as  all  other  LAMBL-type  programs,  also 
creates  graphics  files  called  GRAFXXXX.) 

B.3  The  LAMBL3  program,  which  is  identical  to  LAMBL2,  except  that  it  prints 
only  a  summary  of  the  profile  properties  (item  3  of  the  example  under 
B.4). 

B.4  Sample  run  of  the  LAMBL2  program. 

B.5  Format  of  the  GRAFXXXX  output  files  (these  are  used  for  producing  data 
plots  of  profile  properties). 

B.6  Format  of  the  TOTALXX  output  files  (these  are  used  for  producing  data 
plots  along  x). 


B.l  The  PITOXXXX  Data  File  (Example) 


10  +1 
15  +600 
20  +4 
25  +100 
30  +2.1278 
35  +.1046 
40  +0 

45  +.0001875 
50  +15.45 
55  +.01016 
60  +67 
65  +29 
70  +.008 

120  +1304 

121  +4 

130  +307, +142, 
140  +307, +143, 
150  +310, +143, 
160  +323, +143, 
170  +336, +145, 
180  +351, +150, 
190  +364, +158, 
200  +378,  167, 
210  +390, +178, 
220  +402, +189, 
230  +41 7, +207, 
240  +429, +224, 
250  +444, +248, 
260  +4 56, +272, 
270  +469, +303, 
280  +484, +340, 
290  +495, +389, 
300  +509, $443, 
310  +522, +505, 
320  +536, +568, 


Surface  Code  (eee  Table  I) 

(mm  Hg  abs) 

X  (cm) 

Tq  (deg  F) 

Intercept  of  transducer  calibration  (mm  Hg) 
Slope  of  transducer  calibration  (mm  Hg/count) 
Intercept  of  y— position  calibration  (cm) 

Slope  of  y-position  calibration  (cm/ count) 
Local  surface  pressure  (mm  Hg) 

Probe  diameter  (cm) 

dumber  of  points  in  profile 

Number  of  points  at  the  boundary  layer  edge 

Offset  (cm) 

Profile  Code  (Table  I) 

Number  of  anomalous  points  near  vail 
Beginning  of  data  (y(counts),  p^(counts)) 
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B.2  the  LAMBL2  PROGRAM 
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41 

42 

43 

44 
43 

46 

47 
40 
49 

30 

31 
52 

33 

34 
33 

36 

37 
30 
39 
60 
61 
62 

63 

64 
63 
66 
67 


•10963  •0606434  10.4601  •997765 

•111703  ^0703030  10.7436  *99797 

•11433  .0731107  11.1400  .990277 

•11733  .0760928  11.3953  .998401 

•119768  .0703176  11.9640  .990684 

.122767  .0813037  12.4198  .998786 

•12483  .0035373  12.7327  .990808 

.127433  .0861721  13.1312  .998989 

•129893  .0886003  13.3012  .990409 

•132318  .0912167  13.899Q  .994192 

.134393  .0930065  14.1940  .999293 

•136643  .0953302  14.3267  .499192 

.139268  .0979483  14.9257  .999394 

•  141510  .10«>194  15.2678  .999496 

•  14433  .103002  15.6957  .999^:47 

•  146955  .105623  16.0^51  . 99'»597 

•14883  .107495  16.3805  .999690 

•15103  .110492  16.0371  .999698 

.133705  .147365  17.1223  .999790 

•136703  .115363  17.5794  ,999790 

. i5933  .117406  17.9791  .999790 

•161393  .120047  18.2932  .999094 

. 1 64205  . 1 22859  18.7216  . 999899 

•166643  .125295  19,0929  ,999099 

.16900  .127732  19.4643  I 

.170955  .129608  19,75  1,0001 

.17433  .132983  20.2644  I 


DO  YOU  WISH  A  FILE  HADE  FOR  GRAPHICS  (CR  FOR  YES. I  FOR 


NO>  ?  I 


B.5  The  GRAFXXXX  Files 


These  files  are  outputs  of  the  LAMBL-type  programs  and  contain  the  same 
outputs  as  in  the  printouts  (e.g.  see  B.4)  in  a  form  suitable  for  plotting  on 
a  graphics  plotter.  For  each  point  of  the  profile  given  by  the  XXXX  code,  the 
following  properties  are  provided: 


In  Order 


3uantit> 


N,  i.e.  the  sequential  point  number 


y,  cm 


B.6  The  TOTALXX  Files 

These  files  were  prepared  for  the  express  purpose  of  preparing  graphics 
plots  of  each  integral  or  general  profile  property  vs.  the  streamwise 
distance.  The  XX  designates  the  first  two  symbols  of  the  profile  code  (Table 
II);  e.g.  T0TAL12  would  be  called  to  plot,  say,  0  vs.  x  for  the  smooth  wall 
and  Pq  =  475  mm.  The  following  properties  are  listed  in  these  files: 


In  Order 


Propert\ 


5  ( cm) 

8  (cm) 


Re.S  (cm) 


APPENDIX  C:  DATA  REDUCTION  PROCEDURES 


1 .  Digitization 

Referring  to  Figure  C.l,  the  10-second  data  bursts  (groups)  from  the 
analog  tapes  were  played  back  at  3  3/4  ips  into  the  Rockland  512  Fast  Fourier 
Transform  Computer.  In  emerging  from  the  latter,  the  data  group  was  stored  in 
the  memory  of  the  IBM  9000  Computer  in  digital  form.  Each  such  group  file, 
called  FABCD.BAS,  is  displayed  and  explained  in  Table  C-I  and  is,  of  course, 
the  spectrum  of  the  AC  hot-wire  output  consisting  of  402  numbers  grouped  into 
201  pairs  (intensity  vs.  frequency).  The  upper  frequency  was  320  KHZ  (10  KHZ 
in  playback  time)  and  the  frequencies  contained  in  FABCD.BAS  are  given  as  0, 
1,  2....  (to  be  multiplied  later  by  1.6  KHZ,  which  was  the  real-time  window  of 
the  spectrum).  The  A,  B,  C,  D  in  FABCD.BAS  are  explained  in  the  Table  C-II 
(this  coding  is  the  same  as,  and  an  extension  of,  the  coding  for  the  mean  flow 
data  shown  on  Table  II). 

Once  the  individual  groups  were  digitized,  the  next  task  was  to  collect 
the  groups  belonging  to  each  set  and  compose  a  “masterf ile"  for  that  set. 
This  was  done  with  the  FILEFORM.BAS  program  of  Table  C-III.  These  masterf  iles, 
designated  FABZZ.BAS,  were  used  directly  as  inputs  to  the  data-reduction 
program  STABLE02.BAS  (see  following  section).  Note  the  format  of  these  files 
from  Table  C-I  and  Table  C-IV.  In  the  latter,  a  masterf ile  structure  example 
(for  file  F51ZZ.BAS)  shows  how  hot-wire  "noise”  spectra  (groups)  were 
included.  In  this  example  there  are  two  noise  groups,  F51922  and  F51932  (the 
last  digit  of  these  two  designations  is  not  important  to  this  discussion).  It 
was  found  that  the  data  obtained  for  set  51  (smooth  wall,  free  stream,  Pq  = 
350  torr)  from  x  =  1.3"  to  3.7"  were  subject  to  a  different  electronic  noise 
(in  this  case  F51922)  than  were  the  data  obtained  from  x  =  3.8"  to  6.6"  (this 
happens  because  of  time-dependence  in  the  lab's  electronic  noise  environment 
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and  is  a  rather  minor  point).  Thus  in  this  example  the  net  signal  for  groups 
5113  through  5137  is  obtained  by  appropriately  subtracting  from  them  group 
51922;  these  groups  form  one  "noise  family".  This  is  important  since  the  data 
reduction  program,  when  commanded  to  process  set  F51ZZ,  querries  the  user  on 
the  number  of  such  "noise  families"  (see  following  section).  Table  C-V,  which 
summarizes  this  information  for  all  data  sets,  is  used  as  an  important 
advisory  to  the  computer  operator  during  data  reduction. 

2.  Description  of  the  STABLE02  PROGRAM 

This  program,  written  in  BASIC  language  for  the  IBM  System  9000 
Laboratory  Computer,  is  listed  in  Table  C-VI.  Relevant  information  regarding 
STABLE02  also  appears  in  Table  C-VII  and  Figure  C.2.  The  inputs  to  this 
program  are  the  F — ZZ.BAS  data  files  (Table  C-IV),  which  must  be  in  the  system 
storage  during  the  program  run. 

After  the  LOAD/RUN  command,  the  system  asks  the  operator  seven  questions 
which  are  answerable  by  the  last  seven  numbers  in  each  row  of  Table  C-V  in 
that  sequence.  Note  that  in  the  lines  5-711  the  program  is  set  to  recognize 
the  present  data,  i.e.  modifications  are  necessary  in  order  to  process  future 
data,  say,  at  a  P^  not  used  in  this  test.  Also,  on  three  occasions  (lines 
546,  556,  566),  the  program  splices  together  data  from  two  different  hot¬ 
wires  . 

The  program  will  next  provide  the  option  menu  of  lines  722-1100.  As  soon 
as  the  operator  chooses  an  option  the  program  jumps  to  the  program  line 
indicated  by  lines  1110-1220.  The  function  of  each  option  is  as  follows: 

2.1.  Option  1 :  Raw  Spectrum 

This  is  a  direct  display  of  the  data  as  digitized  by  the  FFT  Computer. 


The  results  are  tabulated  on  the  CRT  screen  only  (volatile). 


2.2.  Opt  ion  2 :  Noise-Free  Spectrum 

This  is  just  like  Option  1,  but  the  resulting  spectral  density  (the  A(f)) 
is  noise-free. 

2.3.  Option  3 ;  Raw  Spectrum  Plot 

Identical  to  Option  1,  but  the  computer  CRT  screen  presents  a  plot  of 
A(f)  vs.  f  (volatile).  Control  of  maximum  f  and  A(f)  in  this  plot  allows  the 
operator  to  magnify  portions  of  the  spectrum. 

2.4.  Opt  ion  4;  Noise-Free  Spectrum  Plot 

As  in  Option  2,  but  a  volatile  net  spectrum  plot  is  obtained  as  opposed 
to  a  tabulation.  In  addition  one  can  make,  on  command,  a  permanent  graphics 
file  of  the  net  spectrum.  All  "graphics"  files  mentioned  herein  are  suitable 
for  making  formal,  publication-grade  plots  by  using  the  Hewlett-Packard  7470A 
System;  they  are,  of  course,  suitable  for  any  additional  analysis  of  the  data 
outside  of  the  STABLE02  program. 

The  data  on  this  separate  file  prepared  in  this  option  are: 

1)  Filename:  SPABCD.BAS 

2)  Coding:  AB:  set  number 

CD:  X  position  in  tenths  of  inch  from  L.E. 

3)  Number  of  variables:  6,  as  follows: 


Computer 


Variable 

Symbol 

Explanation 

Units 

Ist 

N 

Frequency  index  (1-201) 

— 

2nd 

f 

Frequency 

KHZ 

3rd 

F 

Non-dim.  frequency 

— 

4th 

SA 

Raw  spectral  density  (nns) 

volt/1.6  KHZ 

5  th 

NA 

Noise  spectral  density  (rms) 

volt/1.6  KHZ 

6th 

NET 

Net  spectral  density  A(f)  (rms) 

volt/1.6  KHZ 

4)  Number  of  points:  201 


2.5.  Qpt ion  5;  RMS  Variation  Along  Plate 

The  wideband  rms  hot-wire  AC  signal  is  computed  in  this  option  vs. 
position  along  the  plate.  The  results  are  given  by  a  volatile  CRT  screen 
tabulation  or,  on  command,  on  hard  copy  produced  by  the  IBM  9000  printer. 

It  has  been  observed  that  in  the  free  stream  the  wideband  signal  is  very 
close  to  the  wideband  noise.  As  is  frequently  the  case,  the  electronic  noise 
in  the  laboratory  environment  underwent  small  variations  from  one  hour  or  one 
day  to  the  next,  especially  at  the  high  frequencies.  Occasionally  these  small 
variations  made  the  high-frequency  portion  of  the  noise  appear  bigger  than  the 
portion  of  the  signal  spectrum  itself,  producing  a  negative  net  signal  at  the 
high  frequencies.  To  avoid  this,  advantage  was  taken  of  the  fact  that  the 
free  stream  signal  in  the  wind-tunnel  was  insignificant  above  about  100  KHZ. 
Thus,  for  the  "free  stream"  sets  51....  and  61....  the  wideband  rms  signal  and 
noise  were  obtained  by  summing  the  Fourier  components  only  up  to  112  KHZ  (line 
5512  in  the  program). 

Note  that  to  subtract  the  noise  both  for  narrowband  and  wideband  options 
throughout  the  program,  the  usual  procedure  was  obtained  by  assuming  that  the 
noise  and  signal  for  each  Fourier  component  are  uncorrelated  (subtraction  in 
the  squares).  Also,  wideband  rms  signals  were  obtained  by  computing  the 
square  roots  of  the  sums  of  the  squares  of  the  narrowband  Fourier  amplitudes. 

2.6.  Option  6:  RMS  Variation  Along  Plate 

The  only  difference  between  this  and  Option  5  is  the  form  of  the  outputs 
supplied.  Option  6  supplies  a  volatile  graph  of  A(f)  vs.  x  on  the  CRT  screen, 
and  then  asks  if  a  permanent  graphics  datafile  should  be  made.  The  data  for 
such  a  file  are: 


1)  Filename:  RMSAB.BAS 

2)  Coding:  Afi:  set  number 


3)  Number  of  variables:  7,  as  follows: 


Computer 

Variable 

Explanation 

Units 

1st 

X 

Distance  from  L.E. 

inch 

2nd 

X 

Distance  from  L.E. 

cm 

3rd 

REX 

Rex 

— 

4  th 

SQR(REX} 

R  = 

(Rejj)^/2 

— 

5  th 

RETHNOM 

Reg 

nominal  (see  below) 

— 

6th 

RETHACT 

Reg 

actual  (see  below) 

— 

7th 

RMS  I 

Net 

rms 

volts 

4) 

Number  of  Points: 

These  vary  according  to  the  following  scheduli 

SURFACE 

Pq(TORR)  BL/FS 

CODE 

FILENAME  NO. 

OF  POINTS 

SMOOTH 

350  BL 

31 

RMS31 

63 

475 

32 

RMS32 

55 

600 

33 

RMS33 

44 

ROUGH 

350 

21 

RMS21 

72 

475 

22 

RMS22 

65 

600 

23 

RMS  23 

52 

SMOOTH 

350  FS 

51 

RMS  51 

54 

475 

52 

RMS52 

45 

600 

53 

RMS  53 

38 

ROUGH 

350 

61 

RMS61 

65 

475 

62 

RMS62 

29 

The  "NCg  actual"  and  "Re^  nominal"  are  explained  in  the  text. 


2.7.  Option  7 :  Amplitude  Change  With  x  (Tabulations  Only) 

In  this  option  the  net  spectral  density  A(f)  for  a  given  frequency  f  is 
found  as  a  function  of  distance  along  the  plate.  The  results  are  given  both 
in  volatile  form  of  a  CRT  screen  tabulation  and,  if  desired,  also  as  hard 
copy.  (Note:  In  this  option,  as  well  as  in  Option  8,  the  resulting  A(x;f) 
variations  denote  the  directly-measured  experimental  points;  while  in  Option  9 
"theoretical"  values  of  A(f)  additionally  found  by  curve-fitting  the  A(x;f) 


variations  are  also  given.) 


2.8.  Option  8:  Amplitude  Change  With  x  (Graphics) 

This  is  identical  to  Option  7  except  that  the  results  are  displayed  as  a 
volatile  CRT  screen  plot  or,  on  command,  stored  into  a  graphics  file  with  the 
following  characteristics: 

1)  Filename:  AABCDE.BAS 

2)  Coding:  AB:  set  numbers 

CDE:  frequency-index  number  Jl  (1-201,  where  the  actual 
frequency  in  real  time  KHZ  is  given  by  1.6(J1-1)) 

3)  Number  of  points  as  follows: 


Set  No. 

*  No.  of  Points 

Set  No . 

No.  of  Points 

31 

63 

51 

54 

32 

55 

52 

45 

33 

44 

53 

38 

21 

72 

61 

65 

22 

65 

62 

29 

23 

52 

4) 

Number  of  variables:  6 

,  as  follows: 

NO. 

SYMBOL 

EXPLANATION 

UNITS 

Ist 

X 

Distance  from  L.E. 

cm 

2nd 

REX 

Re^ 

— 

3rd 

SQR(REX) 

R 

— 

4  th 

RETHNOM 

Re^  nominal 

— 

5th 

RETHACT 

Reg  actual 

— 

6th 

NET 

A(f)  (spectral  density) 

volts  rms/1.6  KHZ 

2.9.  0] 

Btion  9:  Curve-Fitted  Ami 

ilitude  Chanee  With  x 

The  purpose  of  this  option  is  to  present  the  A(f)  vs.  x  variation  (for 
chosen  f)  of  Options  7  and  8  in  analytic  form,  by  providing  least-squares 
curve-fits  to  the  A(f)  vs.  x  data.  Therefore,  this  option  is  a  "dressing 
room"  in  which  the  degree  of  the  curve-fitting  polynomial  is  chosen  and  a  CRT 
*For  example,  set  31  corresponds  to  file  F31ZZ.BAS. 
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screen  plot  is  quickly  produced  both  of  the  actual  data  points  and  the  fitting 
polynomial.  If  the  fit,  in  the  judgment  of  the  operator,  appears  unsuitable 
then  another  degree  is  chosen  and  the  process  speedily  repeated  until  the  best 
fit  is  achieved. 

This  is  a  very  important  option  of  the  STABLE02  program  where  the  data 
stretch  through  the  laminar  and  transitional  ranges  and  where  preconceptions 
about  the  required  polynomial  fit  can  result  in  significant  errors  in  the 
amplification  rates.  A  more  detailed  account  of  how  this  option  was  used  for 
the  present  data  appears  in  Appendix  D.  This  appendix  notes  that  the  majority 
of  data  were  reduced  using  a  7th-degree  polynomial  fitting  only  the  points  x  < 
’0  cm.  The  latter  range,  of  x,  is  also  optional  during  execution  of  Option  9; 
some  curve-fits  were  made  by  taking  the  entire  range  of  x  (these  are  the  so- 
called  EX - files,  see  below). 

The  option  querries  the  operator  for  the  desired  frequency  to  be 
examined,  the  maximum  x  as  per  above,  the  maximum  value  of  spectral  density 
(amplitude  or  spectral  density  A(f)  in  each  1.6  KHZ  window)  and  the  desired 
degree  of  polynomial  fit  (1-12).  When  the  answers  are  provided,  the  CRT 
screen  shows  a  graph  of  the  data  points  A(f)  vs.  x  for  the  chosen  frequency, 
with  X  ranging  from  0  to  the  limit  given,  and  A(f)  from  0  to  the  amplitude 
limit  given.  The  computer  then  also  gives,  on  the  same  graph,  a  curve  which 
is  the  required  curve-fit.  If  the  operator  wishes  a  different  polynomial 
degree  he  so  commands  until  the  desired  fit  is  reached. 

The  program  next  asks  the  operator  whether  graphics  data  files  of  the 
displayed  result  are  desired.  If  so,  the  operator  is  next  asked  whether  files 
of  (a)  data  points  only,  (b)  amplitude  and  dimensional  amplification,  are 
desired.  The  former  are  designated  by  ”PT"  and  the  latter  by  "CF".  The  PT 
files  are  set  up  as  follows: 


1)  Filename:  PTABCDE.BAS 


2)  Coding:  AB:  set  number 

CDE:  frequency  index  (1-201).  Actual  frequency  is 
f(KHZ)-1.6  CDE. 


3)  Number  of  points:  Depends  wholly  on  the  desired  range  of  x,  as 
stated  during  execution  (program  line  9530).  Since  the  CRT  display  of  this 
file  takes  place  before  a  decision  to  make  such  a  file  is  taken,  the  operator 
can  count  the  data  points  included  in  the  graph.  In  analyzing  the  present 
test  results,  the  x  range  was  limited  to  10  cm  for  curve-fitting  purposes;  in 


0  <  X  <  10  cm 

the  points  in 

FT  f  iles  were: 

Set  No. 

No.  of  Points  Set  No. 

No. 

of  Points 

31 

36 

51 

27 

32 

37 

52 

31 

33 

36 

53 

30 

21 

36 

61 

29 

22 

33 

62 

29 

23 

37 

4)  Number  of  variables:  3,  as  follows: 

Variables 

Symbol 

Exolanation 

Units 

let 

XCF 

Distance  from  L.E. 

cm 

2nd 

RETHNOM 

Re  ^  nominal 

— 

3rd 

YCF 

Spectral  density  (amplitude  A(f)) 

vrms/1.6  KHZ 

Recall 

that  the  A - 

-  files  of  Option  8  produce  much 

the 

same 

type  file 

with  the  total  number  of  x  stations  recorded  for  each  set. 

The  files  CF  produce  curve-fits  of  the  data  displayed  on  the  CRT  screen 
during  this  option  and  stored  in  the  FT  files.  All  data  contained  in  this 
file  are  results  of  the  curve  fit;  for  example  the  3rd  variable  in  the  FT  file 

(see  above)  is  the  experimentally  determined  A(f),  whereas  CF - will  produce 

a  A(f)  which  is  the  theoretical  A(f)  lying  on  the  polynomial  curve.  In 
addition  CF  produces  the  dimensional  amplification  rate  dA/dx  as  well  as  the 
usual  non-dimensional  one: 


2A 
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1)  Filename:  CFAfiCDEF.BAS 

2}  Coding:  AB:  Set  number 

CDE:  frequency  index  (1~201) 
F:  polynomial  degree 

3)  Number  of  points:  100 

4)  Number  of  variables:  5,  as  follows: 


Units 

cm 


1st  G(Nl)  Distance  from  L.E.  cm 

2nd  R  R  = 

3rd  P{N1)  Spectral  density  (theoretical  A(f))  vrms/1.6  KHZ 

4th  Pl(Nl)  Dimensional  amplification  rate  vrms/cm 

5th  —  Non-dimensional  amplification  rate  - 

NOTE:  In  analyzing  the  present  data,  this  option  was  occasionally  also  run 
using  the  full  available  range  of  x  (this  was  suitable  for  studying 
amplification  rates  of  the  second  mode).  In  this  case,  the  polynomial  degree 
was  again  7  and  the  file  was  prefixed  "EX”  rather  than  "CF";  otherwise 
everything  was  identical  to  "CF"  files.  In  this  case,  if  one  wished  to  plot 
both  the  A(f)  points  and  their  curve-fits  on  the  same  graph,  one  produced  and 

used  A - files  from  Option  8  together  with  EX - files  (just  as  one 

matched  PT -  and  CF -  files  for  x  <  10  cm). 

2.10.  Option  10 :  Poles  of  the  Stability  Diagram 

The  purpose  of  this  option  is  to  provide  points  in  the  stability  diagram 
marking  the  "lower"  and  "upper"  neutral  branches  as  well  as  the  maxima  and 
minima  in  the  amplification  rates.  This  is  done  exclusively  using  the  curve- 
fitted  curves  of  -«^vs.  x  (or  R,  etc.).  The  lower  branch  is  defined  as 

the  point  on  the  F,  Re^  plane  where  -otXat  a  fixed  frequency)  changes  from 
negative  to  positive,  and  the  opposite  is  true  for  the  "upper  branch." 

The  results  are  presented  as  a  volatile  CRT  screen  tabulation,  after 


^>1  ■.»  \\  r j  v  rj>  rr  v  v  ■v.-’.n- 


which  the  operator  can  obtain  hard  copy  of  the  tabulation  or  ask  the  machine 
to  prepare  any  of  four  graphics  files  as  follows: 

FILENAME  PURPOSE 

LOWERAB.BAS  Locates  the  lower-branch  points  in  the  diagram. 

UPPERAB.BAS  Locates  the  upper-branch  points  in  the  diagram. 

MAEAMAB.BAS  Locates  maxima  in  the  -  otj^vs.  R  curve. 

MINAMAB.BAS  Locates  minima  in  the  -ot^vs.  R  curve. 

AB:  set  number 

Thus  each  data  set  produces  four  files,  as  above.  Also  note  that  for 
each  of  these  files,  at  each  frequency  there  can  be  more  than  one  point  in 
each  file,  e.g.  there  may  be  two  or  more  upper  neutral  branches  (this  was,  of 
course,  a  key  issue  in  this  experiment  where  the  data  disclosed  the  existence 
of  more  than  one  unstable  region). 


The  number  of  points  were  as  follows: 


Set  No.  (AB) 

LOWERAB 

UPPERAB 

MAXAMAB 

MINAMAB 

31 

184 

113 

198 

271 

32 

161 

67 

127 

217 

33 

129 

73 

129 

180 

21 

224 

153 

203 

266 

22 

159 

110 

192 

207 

23 

112 

88 

200 

198 

51 

303 

281 

288 

317 

52 

234 

236 

279 

275 

53 

199 

189 

261 

270 

61 

375 

284 

291 

381 

62 

379 

290 

294 

385 

VARIABLES : 

6 

6 

7 

7 

As  indicated,  the  number  of  variables  differ,  and  were  as  follows: 


HI 


Variable 


1st  flKHZ;  Frequency  KHZ 

2nd  F  Non-dimensional  frequency  - 

3rd  XL(J)  Distance  from  L.E.  cm 

4th  SQR(REX)  R 

5th  RETHNOM  Rcg  nominal 

6th  RETRACT  Re^  actual 

and  for  MINAMAB,  MAXAMAB  only, 

7th  PIMAX  (-a.)  for  MAXAM 

A  mAxxinuni 

7th  PIMIN  (-«.)  for  MINAM 

•i  tDinunusi 

2.11.  Option  11:  Amplification-Rate  Spectrum 

The  purpose  of  this  option  is  to  produce  amplification  spectra,  i.e.  the 
variation  -or-vs.  F  for  chosen  values  of  R  (or  x  or  Re  ).  Restrictions  are: 
(a)  the  option  is  based  on  7th  degree  polynomials,  (b)  the  7th-degree  fit  is 
again  done  only  for  x  <  10  cm,  (c)  the  amplification  spectra  for  R  higher  than 
those  corresponding  to  x  »  10  cm  will  therefore  be  possible  but  unreliable, 
(d)  the  spectra  are  limited  below  F  corresponding  to  160  KHZ. 

The  option  produces  volatile  CRT  screen  tabulation  of-oi^vs.  F.  On 
command,  graphics  files  are  produced  as  follows: 

1)  Filename:  ASABCDE.BAS 

2)  Coding:  AB:  set  number 

CDE:  R  value 

3)  Number  of  points:  100 

4)  Number  of  variables:  2,  as  follows: 


Variable 


Explanation 


F  Non-dimensional  frequency  F 

RATE  (Jl)  Non-dimensional  amplification  rate  -<x. 


Option  12:  Total  Amplification  Spectrum 


The  purpose  of  this  option  is  to  find  the  total  amplification  of  the 
disturbance  of  a  given  frequency  f  between  two  points  and  X2  on  the  plate 
or,  rather,  Rj  and  R2*  The  computer  does  this  with  the  curve-fitted  A(R) 
variations,  first  by  calculating  these  variations  with  Option  9  (again  curve¬ 
fitting  with  7th  degree  polynomials  the  data  in  x  <  10  cm).  The  initial 

R  =  Ri  is  chosen  to  be  150  (this  is  near  the  minimum  x  measured,  near  and 
0  1 

below  which  the  curve-fits  are  unreliable).  Then,  for  an  operator-chosen  R2 
the  machine  provides  A(R2)/A(R]^).  Note  that  R][  can  be  changed  in  line  12720 
of  the  program. 

Graphics  files  are  also  made  as  follows: 

1)  Filename:  TOTABCDE 

2)  Coding:  AB:  set  number 

CDE:  R2  value 

3)  Number  of  points:  100  (to  160  KHZ) 

4)  Number  of  variables:  2,  as  follows: 


Variable 


Explanation 


Non-dimensional  frequency 

A(f;R2)/A(f;Ri) 


Units 


2.13.  Changes  in  the  Curve-Fitted  x-Range 

If  desired,  the  curve-fitted  x  range  can  be  easily  changed.  For  example 
to  curve-fit  from  0  to  15  cm.,  do  the  following: 

For  Option  9:  Answer  querry  (line  9530)  with  "15" 


For  Option  10 
For  Option  11 
For  Option  12 


Change  line  10540  to  read:  LET  Xl  =  15 
Change  line  12110  to  read;  LET  XI  -  15 
Change  line  12630  to  read:  LET  XI  =  15 


TABLE  C-I 


A.  FILES  : 
Numbers  in 


B.  FILES  ■ 
Numbers  in 


FORMAT  OF  F - .BAS  AND  F--ZZ.BAS  FILES 


' - .BAS  (GROUP  DATA) 


Sequence  (Example): 


3214 

Group  Code 

0 

First  frequency 

First 

59.86 

First  A(f) 

pair 

1 

Second  frequency 

Second 

62.59 

Second  A(f) 

pair 

2 

Third  frequency 

201 

pairs 

(402 

numbers) 


403 

numbers 


200  Last  frequency  ~|  Last 

0.135  Last  A(f)  J  pair 


■--ZZ.BAS  (SET  DATA) 


Sequence  (Example) 


3202 

Group  Code  1 

0 

First 

frequency 

53.2 

First 

A(f) 

First 

1 

Second  frequency 

group 

62.12 

Second  A(f)  J 

3203 

Group 

code 

0 

First 

frequency 

Second 

55.7 

First 

A(f) 

group 

3216 

Group 

code 

0 

First 

frequency 

Last 

50.2 

First 

A(f)  i 

group 

3281 

Group 

code 

0 

First 

frequency 

First 

1.32 

First 

A(f) 

noise 

1 

Second  frequency 

group 

0.97 

Second  A(f) 

3282 

Group 

code 

Second 

0 

First 

frequency 

noise 

1 .02 

First 

A(f) 

group 

Set  masterfile 
F32ZZ.BAS  (Number 
of  numbers  depends 
on  number  of  groups 
and  noise  groups.). 


TABLE  C-III.  THE  FILEFORM  PROGRAM 

Program  FILEFORM: 

6  DIM  Y(55,  405), 1(100), J(405) 

10  INPUT  "TYPE  DESIRED  SET  MASTERFILE  NAME  (IN  QUOTES )";A$ 

15  IF  A$="Q"  GOTO  10000 

20  OPEN  A$  FOR  OUTPUT  AS  FILE  #1 

30  FOR  1=1  TO  100 

40  INPUT  "TYPE  INPUT  GROUP  FILENAME  IN  QUOTES  ('END'  IF  DONE)";B$ 
50  IF  B$="END"  GOTO  9000 
60  OPEN  B$  FOR  INPUT  AS  FILE  #2 
70  FOR  J=1  TO  403 
80  INPUT  #2,Y(I,J) 

90  PRINT  #1,Y(I,J) 

100  NEXT  J 
105  CLOSE  n 
110  NEXT  I 
9000  aOSE  #1 


10000  END 


TABLE  C-IV 


DATA-SET  FILE  FORMAT  ACCEPTABLE  TO  THE 
STABLE02  PROGRAM 


51 B,?  =  F51922 
"NOISE  GROUPS" 
5,33  =  F5132 


TOTAL  GROUPS 


(1):  3rd  question  of  program,  (?)  4th  question, 
(3):  5th  question,  (4);  1st  question 


TABLE  C-V 


SUMMARY  OF  SET  DATA  PROCESSED  BY  THE 
DATA-REDUCTION  PROGRAM  STABLE02 


BOUNDARY  LAYER  DATA 

SET  TOTAL  NO.  OF 


LAST  x(i); 

1st  2nd  3rd 


NO.  OF 


SURFACE 

PqCTORR) 

NAME  GROUPS 

FREQUENCIES 

FAMILY 

FAMILY 

FAMILY 

FAMILIES 

SMOOTH 

350 

F31ZZ^2) 

64 

201 

61 

100^3) 

100^3^ 

1 

475 

F32ZZ 

57 

201 

29 

57 

100 

2 

600 

F33ZZ 

45 

201 

47 

100 

100 

1 

ROUGH 

350 

F21ZZ 

75 

201 

43 

72 

77 

3 

475 

F22ZZ 

67 

201 

37 

73 

100 

2 

600 

F23ZZ 

53 

201 

56 

100 

100 

1 

FREE  STREAM  DATA 

SMOOTH 

350 

F51ZZ 

56 

201 

37 

66 

100 

2 

475 

F52ZZ 

47 

201 

37 

53 

100 

2 

600 

F53ZZ 

40 

201 

37 

47 

100 

2 

ROUGH 

350 

F61ZZ 

66 

201 

77 

100 

100 

1 

475 

F62ZZ 

31 

201 

17 

37 

100 

2 

Notes:  1)  Indicates  last  group  in  set,  in  tenths  of  an  inch  from  L.E. 

2)  Terminator  .BAS  needed  to  process  by  STABLE02 

3)  The  "100"  is  a  code  with  no  physical  meaning 
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TABLE  C-VII 


OUTPUTS  PROVIDED  BY  STABLE02.BAS 


OPTION  PURPOSE 


CRT*  CRT* 

TABLE  PLOT  HARDCOPY 


NAHE(S)  OF 
CREATED  FILE(S) 


Raw  apectrum,  each  group  % 

Net  spectrum,  each  group  x 

Raw  spectrum,  each  group 

Net  spectrum,  each  group 

RMS  variation  with  x,  x 

entire  set 

RMS  variation  with  x, 
entire  set 

Amplitude  variation  with  x,  x 

single  frequency 

Amplitude  variation  with  x, 
single  frequency 

Curve-fits  of  single-frequency 
amplitude  variation  with  x 


SP - .BAS 


RMS— .BAS 


10  Stability  Diagram  Poles 


LOWER— .BAS 
UPPER— .BAS 
MAXAM— .BAS 
MINAM— .BAS 


11  Amplification  rate  spectrum 
at  some  fixed  x  (or  Re,^) 

12  Total  amplification  between 
Rg*  150  and  given  R 


TOT - .BAS 


*Volatile  outputs 


ISO 


3  3/4  IPS  PLAYBACK 


ANALOG  TAPE  RECORDER 
HONEYWELL  7600 
CHANNEL  7  OR  2 


I  COMPILER 
I  PROGRAM 
"FILEFORM" 


1024  AVERAGES 
GAIN  8-32  VOLTS 
SENSITIVITY  0.1-1 
RANGE  10  KHZ  (PLAYBACK) 
320  KHZ  (REAL) 


FILES  FABCD.BAS 
AB:  SET  CODE 

CD:  X  (TENTHS  OF  INCH  FROM  L.E.) 


FILES  FABZZ.BAS 
AB:  SET  CODE 


-  END  OF  DIGITIZATION 
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OVERVIEW  OF  STABLE02  PROGRAM 
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C.2  Overview  of  the  STABLE02  Program. 


1  o<a _ 


The  important  variables  in  the  discussion  of  stability  are  the  net 
spectral  density  e^(f ,R)  and  the  amplification  rate  -  Q(^(f ,R).  The  e^Cf.R) 
is  measured  as  mean-square  volts  (VMS)  within  a  narrow  frequency  window  which 
in  our  case  was  1.6  KHZ.  In  the  present  experiment  too,  the  data  were 

A 

collected  and  are  presented  as  the  square  root  of  e^(f,R),  which  is  called 
A(f)  or  simply  A.  The  spectra  of  the  amplifier  output  at  a  given  R  are  then, 
in  this  test,  given  by  A(f)  vs.  f.  The  variation  A(R)  at  some  fixed  frequency 
f  will  be  called  the  "amplitude  change"  and  can  also  be  alternately  shown  as  A 
vs.  X  or  vs.  Reg  ,  or  versus  any  other  such  variable  connected  with  variations 
along  X. 

In  order  to  illustrate  the  role  of  the  flow  as  an  amplifier,  most 
stability  discussions  and  comparisons  with  theory  deal  with  the  amplification 

-<v, 

^  2a  dR. 

Thus  depends  on  taking  derivatives  of  curves  formed  by  data  points,  which 
introduces  inaccuracies  deserving  a  brief  discussion  of  past  work.  The  data 
of  Demetriades  (References  6,  7  and  8)  and  Stetson  (Reference  9)  were  obtained 
by  first  measuring  the  spectra  A(f)  at  a  series  of  R  values,  crossplotting 
these  vs.  R  at  fixed  f,  and  finally  passing  some  sort  of  continuous  curve  A(R) 
through  the  points  in  the  A-R  plane  from  which  dA/dR  was  computed. 
Demetriades  (References  6  and  7)  and  Stetson  et  al.  (References  9,  10  and  11) 
chose  4th-degree  polynomials  to  curve-fit  their  points  with  the  least-squares 
method.  Laufer  and  Vrebalovich  (Reference  5)  and  Kendall  (Reference  12)  do 


not  specify  the  curve-fitting  method.  In  their  caset  however,  it  should  be 
noted  that  amplitude  data  A(R)  vs.  R  were  available  in  the  form  of  a 
continuous  line,  i.e.  an  "infinite"  number  of  points  in  the  A-R  plane  (this 
was  done  by  plotting  directly  a  single  Fourier  component  output  vs.  distance 
on  an  z-y  plotter).  Such  a  large  number  of  points  should  improve  the  accuracy 
of  measuring  slopes,  although  "scatter"  can  exist  on  a  continuous  line  also. 

In  the  present  work  the  A(R)  data  were  least-squares  curve-fitted  as  in 
References  8  and  9  but  with  considerable  reservations  which  apply  in 
retrospect  to  these  references  as  well.  If  a  set  of  data  points  is  not  too 
badly  scattered  and  defines,  to  the  naked  eye,  a  fairly  simple  curve,  then  it 
is  true  that  just  about  any  higher-degree  polynomial  will  fit  the  set  well. 
Even  then,  however,  the  slope  of  any  fitted  curve  will  depend  on  (a)  the 
polynomial  degree  and  (b)  the  number  of  points  included.  This  dependence  can 
become  intolerable  if  there  is  data  scatter  and/or  the  data  imply  a  "complex" 
curve  (e.g.  several  inflection  points).  Ho  justification  has  been  given  in 
References  8  and  9  as  to  why  the  polynomial  degree  was  chosen  at  4,  and 
therefore,  it  is  not  easy  to  determine  the  accuracy  of  their  results  (Figure  6 
of  Reference  8,  Figures  16  of  Reference  9). 

Even  if  the  curve-fits  for  the  calculations  in  References  8  and  9  were  to 
be  justified  by  the  evidently  "good  behavior"  of  the  points  to  be  fitted, 
curve-fitting  of  the  present  data  looked  much  more  difficult  at  the  beginning. 
For  example,  in  Figure  D.l,  we  see  that  there  are  as  many  as  three  extrema 
defined  by  the  A(R)  (in  this  case  the  A(x))  curves,  even  in  the  limited  range 
0  <  X  <  10  cm.  Furthermore,  in  most  cases  the  large  signal  increases  at  high 
X  tended  to  so  dominate  the  A(R)  picture  that  any  curve  fitting  it  would  fit 
poorly  the  small  but  significant  changes  seen  at  small  x.  It  therefore 
appeared  unreasonable  to  make  a  single  curve-fit  of  the  entire  x  range  or  to 
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favor  some  polynomial  degree  in  advance*  It  was  instead  decided  to  study  in 
some  detail  the  regions  into  which  x  should  be  divided  for  individual  curve- 
fits  and  the  degree  fitting  the  data  best.  To  avoid  biasing,  however,  it  was 
decided  to  use  the  same  polynomial  degree  for  all  curve-fits  once  this  proper 
degree  was  determined.  The  judgment  of  the  latter  was  made  by  plotting  the 
polynomial  slopes. 

The  Effect  of  the  Polynomial  Degree 

The  search  for  the  proper  degree,  shown  by  examples  on  Figures  D.1-D.3, 
was  done  by  first  looking  at  the  data  in  0  <  x  <  10  cm,  i.e.  mostly  in  the 
"first  mode"  region.  In  this  region  the  figures  show  the  data  fitted  by  six 
polynomials  of  degree  2-8.  In  each  example  of  Figures  D.1-D.3  there  are  two 
plots,  the  first  plotting  A  (i.e.  spectral  density)  vs.  x,  and  the  second  the 
amplification  rate  (actually  dA/dx)  vs.  x.  Each  polynomial  was  found  solely 
by  using  the  points  in  0  <  x  <  10  cm,  even  though  one  example  in  Figure  D.3 
presents  a  detail  of  the  curves  confined  into  0  <  x  <  4  cm.  Note  that  the 
examples  are  so  chosen  as  to  include  different  types  of  amplitude  changes 
ranging  from  near-monotonic  amplification  (e.g.  f  =  40  KHZ,  smooth  wall,  p^  * 
475  mm)  to  near-monotonic  damping  (e.g.  f  *  76.8  KHZ,  350  mm). 

Judging  especially  from  the  amplification  rate  plots,  it  is  seen  that  the 
polynomial  degree  has  to  be  rather  high  (6-8)  before  the  curves  begin  to 
settle  toward  a  shape  nearly  independent  of  it.  After  studying  such  curves  it 
was  decided  to  curve-fit  with  7th  degree  polynomials. 

The  Effect  of  the  X-Range 

Once  the  suitable  polynomial  degree  was  "narrowed  down",  a  test  of  the 
appropriate  x-range  was  made  with  plots  such  as  those  on  Figures  0.4  and  D.5. 
The  procedure  here  was  as  follows.  First,  the  polynomial  degree  is  chosen 
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(only  degrees  7  and  8  are  used  since  those  were  the  most  promising;  see 
above).  Then  the  points  in  the  x  range  0-5  cm  are  fitted  with  such  a 
polynomial.  The  procedure  is  then  repeated  for  the  ranges  0-10  and  0-20  cm. 
The  results  are  again  plotted  in  the  A-x  and  dA/dx-x  planes. 

Figure  D.4  shows  a  case  where  in  the  range  1  <  x  <  4  cm  (a  critical  range 
of  the  first-mode  region)  the  amplification  rates  are  not  sensitive  to  the 
polynomial  degree  choice  (7  or  8),  whether  the  points  in  the  range  5  <  x  <  10 
cm  are  included  or  not.  When  the  points  in  10  <  x  <  20  cm  are  included  in  the 
curve-fit,  then  the  curve-fit  seems  improper.  On  Figure  D.5  another  case  is 
shown  where  the  insensitivity  to  the  polynomial  degree  remains  but  where  a 
consideration  only  of  the  points  in  the  range  0  <  x  <  5  cm  produces  an 
inappropriate  curve-fit. 

On  the  basis  of  this  test,  it  was  decided  that  the  curve-fits  needed 
would  consist  of  7th  degree  polynomials  fitting  the  range  0  <  x  <  10  cm.  This 
was  the  method  employed  for  the  majority  of  the  figures  in  the  text  and  for 
the  bulk  of  the  analysis  and  discussion.  Seventh-degree  polynomials  fitting 
the  range  0  <  x  <  20  cm  were  also  used  for  a  small  segment  of  the  text  where 
it  was  desired  to  study  and  discuss  the  second-mode  region. 


spectral  density  CVRMS/1.6KHZ>  XIO'-^-SPECTRAL  density  (VRMS/1.  SKHZ) 
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D.2  Curve- fit  examples  at  57.6  and  76.8 
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D.3  Curve-fit  examples  at  88  and  126,4 
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D.4  Influence  of  included  x  range  on 
curve-fits. 
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D,5  Influence  of  included  x  range  on 
curve-fits . 


